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INTRODUCTION 


The lobe mixer is a device currently being used for thrust augmentation 
on a variety of turbofan engines. The forced mixer is designed to mix 
primary (core) and secondary (fan) flow streams before they enter the exhaust 
nozzle, thus providing the nozzle with air having a more uniform energy which 
results in a more uniform velocity at the nozzle exit. It has long been 
known that by mixing the fan and turbine exhaust streams of a turbofan engine 
in this manner prior to expansion through the exhaust nozzle, a small but 
significant performance gain may be realized. The level of these gains 
depends on trade-offs between the degree of mixing of the two streams and the 
viscous losses incurred in the mixing process. 

To date, the performance of lobe mixers has been determined almost 
entirely through experiments. These experiments are sufficient to determine 
the relative merits of one mixer configuration as compared to another, but 
in the absence of detailed understanding, it is both expensive and time 
consuming to develop unproved lobe mixers in this manner. The overall mixer 
flow field consists of two major flow regions: (1) the flow upstream of the 

lobe exit plane (within the lobes) and (2) the flow downstream of the lobe 
exit plane where the mixing actually takes place. The present study concen- 
trates upon the actual mixing process downstream of the lobe exit plane, 
through development of a computational procedure applicable in this region 
and capable of predicting the three-dimensional mixing process in detail. 

The present analysis and computational procedure represent a capability 
which can provide further understanding of the mechanism by which mixers 
produce a gain in performance. Ultimately, the present analysis may enable 
designers to predict practical mixer flows in detail and to screen or otherwise 
evaluate alternative mixer configurations. 

Although the geometry of the mixing duct Itself is axlsymmetrlc and thus 
effectively contains only radial and axial variations, the flow entering the 
mixing duct is normally three-dimensional and consists of two streams, the 
cold fan stream and the hot core stream. The inlet flow contains free shear 
layers and thermal mixing layers downstream of the mixer lobes. Since the 
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iBlxer lobe geonetry forces en azimuthal variation in the flow, the mixing 
layers of primary importance In the present analysis are Inherently three- 
dimensional, and thus various methods of analysis based on two-dimensional 
theories were excluded from consideration. The flow field downstream of the 
lobe exit plane and within the mixing duct can be treated, however, within 
the general framework recently developed by Briley & McDonald (Ref. 1) for 
the prediction of three-dimensional viscous subsonic flows with large 
secondary velocities. The analysis of Briley & McDonald was developed for 
three-dimensional subsonic viscous flows at high Reynolds number which have 
a predominant primary flow direction with secondary flow In transverse planes. 
A key restriction of the analysis is that once a primary flow direction is 
established, reversed flow in the primary flow direction (flow separation) 
does not occur. The analysis synthesizes concepts from potential flow theory, 
secondary flow theory, and from an extension of three-dimensional boundary 
layer theory in a manner which allows efficient solution by forward-marching 
numerical integration techniques, given only knowledge of the potential flow. 
Repeated and costly (iterative) sweeping of the entire three-dimensional flow 
field is not necessary to obtain a valid approximation to the flow. As a 
consequence, numerical solution of the derived governing equations generally 
requires considerably less computer time than would a solution of the full 
three-dimensional Navler-Stokes equations. An a priori potential flow solution 
nevertheless accounts at least approximately for the important elliptic 
Influence of downstream boundary conditions and determines the primary flow 
direction. In Ref. 1, the analysis was employed to compute laminar flow in 
curved rectangular ducts and a turbine blade passage at high Reynolds number 
but relatively low ^!ach number. In that application, the strong secondary 
flows and comer vortices associated with the endwall regions and caused by 
large turning angles of the passage, were of primary interest. The present 
application to lobe mixer flows Includes a two-equation turbulence model, 
and solves an energy equation to account for thermal effects arising from 
large temperature differences between the hot turbine exhaust and cool fan 
streams. The analysis is designed for Internal flows which have a Mavlnwnn 
average Mach number less than 1.0 at each axial station, although Mach numbers 
exceeding unity may occur locally within the flow field. 
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ANALYSIS 


General Approach 

The present analysis includes a development in a general vector form 
which follows Ref. and other material which pertains to the application 
to turbulent flow in lobe mixers. The analysis is based on approximations 
made relative to a curvilinear but not necessarily orthogonal coordinate 
system fitted to and aligned with the flow geometry under consideration 
(cf. Fig. 1). The coordinate system is chosen such that the streamwise or 
marching coordinate either coincides with or is at least approximately 
aligned with a known inviscid primary flow direction, as determined for 
example by a potential flow for the given geometry. Transverse coordinate 
surfaces must be either perpendicular or nearly perpendicular to solid walls 
or bounding surfaces, since diffusion is permitted only in these transverse 
coordinate surfaces. One obvious coordinate system would take the local 
direction of the velocity from an inviscid potential flow as the marching 
coordinate, and curvilinear surfaces of constant potential as transverse 
coordinate surfaces, although such a "potential flow coordinate system" is 
not essential (and may not even be convenient) . 

Equations governing streamwise vorticlty and a scalar viscous correction 
u^ to a known inviscid primary flow velocity are derived utilizing assump- 
tions which permit forward-marching solution, provided reversal of the 
composite streamwise velocity does not occur. Terms representing diffusion 
normal to transverse coordinate surfaces are neglected. Approximate pressure 
gradients are derived from the inviscid primary flow and Imposed in the 
streamwise momentum equation. These pressure gradients are the sole means of 
accounting for the elliptic influence of downstream boundary conditions in 
curved flow geometries, although they may be corrected as part of the forward- 
marching process if a suitable correction is available. Secondary flow 
velocities are determined by scalar and vector potentials in the transverse 
coordinate surfaces, to suppress streamwise elliptic behavior requiring 
downstream boundary conditions. 

The use of streamwise vorticlty in obtaining secondary flow velocities 
avoids the use of transverse momentum equations and hence the major Influence 
of transverse pressure gradients in {Sredictlng secondary flow velocity. 
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As a consequence, the method does not require or yield a direct prediction 
for the overall viscous pressure field. The viscous pressure field can never- 
theless be determined ^ posteriori using the computed composite velocity as 
input for a scalar three-dimensional solution of the Poisson aquation for 
pressure, obtained as the divergence of the vector momentum equation. Although 
the presence of reversed flov in the primary flow coordinate direction (flow 
separation) normally precludes the use of forward-marching solution procedures 
as described here, special treatments which exclude the separated region from 
the computational domain are both possible and useful in some applications. 

This general approach has been applied to flow in smoothly curved internal 
flow passages with large turning by Briley & McDonald [1 & 2], to the 
airfoil tip vortex generation process by Shamroth & Briley L3]> and is 
also applicable to the forced lobe mixer problem considered here. 


Primary-Secondary Velocity Decomposition 

The analysis is based on decomposition of the overall velocity vector 
field U into a primary flow velocity and a secondary flow velocity U^. 
The overall or composite velocity is determined from the super-position 


U - Up + u. 


The primary flow velocity is represented as 


Up - Uju, 


( 2 . 1 ) 


( 2 . 2 ) 


where is a known Inviscld primary flow velocity satisfying slip conditions 
and determined for example from an a^ priori potential flow solution in the 
geometry under consideration. The (non-dimensional) scalar quantity u^ is a 
viscous velocity profile factor which introduces viscous shear layers and may 
also correct for internal flow blockage effects. In the present case of the 
lobe mixer u^ accounts for both boundary layers and the free shear layer which 
exists between the fan and turbine streams. The viscous velocity correction 
u^ is determined from solution of a primary flow momentum equation. As will 
become apparent, for the lobe mixer it is convenient to solve for the 
component of D^ in the primary flow direction rather than u^. 
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To suppress^ 8 treamwlse elliptic behavior, the secondary flow velocity 
is defined by and thus presumed derivable from scalar and vector surface 
potentials denoted ^ and respectively. If i^ denotes the unit vector 
normal to transverse coordinate surfaces (also presumed here to be in the 
direction of the marching coordinate) , if p is density, and if is an 
arbitrary constant reference density, then is defined by 

where is the surface gradient operator defined by 

It follows that 1^ • 0 and thus lies entirely within transverse 

coordinate surfaces. Equation (2.3) is a general form permitting both 
rotational and irrotatlonal secondary flows and will lead to governing equa< 
tlons which are elliptic in transverse coordinate surfaces and which are 
therefore solvable within a forward-marching context. The overall velocity 
decomposition (2.1) can be written 

U ■ U u + 7j<^ + ( Pq //>)7xt 

(2 

“ + ^Po ^P^ ^ ^ 


(2.3) 


(2.4) 


Surface Potential Equations 

Equations relating ^tid ijj with u^, p, and the streamwlse vortlcity 
component 12^ can be derived from £q. (2.5) as follows: 


since the last term in each of Eqs. (2.6, 2.7) is zero by vector identity, 

Eqs. (2.6, 2.7) can be written as 

• -V-pUjUy (2.8) 

(2.9) 

which are elliptic for i and 4* in transverse coordinate surfaces. The last 

term in Eq. (2.9) is Identically zero in a potential flow coordinate system 

for which 1^ and have the same direction, and is small if i^ and are 

approximately aligned. Given a knowledge of u^, G, and p, the surface 

potentials ^ and i|) can be determined by a two<-dimensional elliptic calculation 

in transverse coordinate surfaces at each streamwlse location. In turn. U 

^ s 

can be computed from Eq. (2.3), and the composite velocity U will satisfy 
continuity. Equations for u^ and are obtained from the equations governing 
momentum and vortlclty, respectively. 

Primary Momentum and Pressure Approximation 

The streamwlse momentum equation is given by 

?, •[(U*V)U +(Vp)//>] - I,* F (2.10) 

where p is pressure and pF is force due to viscous stress. Terms in F 
representing streamwlse diffusion are neglected; however, since the viscous 
terms are complex for compressible flow, the modified viscous force is 
temporarily denoted F', and further consideration of viscous terms is deferred 
to a later section. 

The remaining assumption for Eq. (2.10) concerns the pressure gradient 
term and is designed to permit numerical solution as an initial value problem. 
An obvious pressure approximation for curved flow geometries is to impose 
pressure gradients from an inviscld potential flow (Briley [4]), and for 
internal flows, to correct these with a mean pressure gradient term depending 
only on the x^ coordinate (Patankar t Spalding [5]). This approximation is 
appropriate both for flows consisting of an irro rational core region with 
thin shear layers on bounding surfaces, and also for some fully viscous flows. 
Depending on the particular application, approximate downstream boundary 
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conditions such as the commonly used "Kutta condition" may be required, and 
periodicity conditions may be imposed to simulate cascades of airfoils. This 
approximation invokes the boundary layer concept that, to first order, shear 
layers do not alter inviscid pressure gradients. 

A somewhat more general approximation (in that need not be a potential 
flow) can be derived, and the consequences of such an assumption are clarified 
on examination of the inviscid momentum equation written as 

Vp + /)V(q^/2) ■ /)U X (2.11) 

where ■ U • U, and vorticity fi ® V x U. In an orthogonal coordinate system, 
the i^ component of Eq. (2.11) is 


l,*[vp + /?V(q2/2)] - />(u2il3-U3fl2) 


( 2 . 12 ) 


For small transverse vorticity £2^ and small transverse velocities u^, 
(relative to the coordinate system being used) , the right-hand side of 
Eq. (2.12) is of second order and may be neglected. This small shear assump- 
tion is familiar as part of the "secondary flow approximation" (cf. Hawthorne 
[6]), although it should be emphasized that the present analysis accounts for 
"distortion of Bernoulli surfaces" by means of u^. 

If the right side of Eq. (2.12) is neglected, if pressure gradients are 
derived from the inviscid velocity U^. by setting q^ * • Uj, and if a mean 

viscous pressure correction p^(x^) is introduced, the pressure approximation 
can be written as 


i,*Vp * i,-[vp^(x,)- />7(Uj-Uj)/2] 


(2.13) 


where p^ is the Imposed pressure. Typically, for Internal flows, p^ is 
determined to ensure that an Integral mass flux condition Is satisfied, such 


as 


J I, • />UdA 


constant 


(2.14) 
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For external flows, ■> 0. If an Internal flow has an average Mach number 
of 1.0 and thus becomes choked, it Is not possible to determine a value for 
p^ such that Eq. (2.14) is satisfied. In this Instance, the solution down- 
stream is indeterminant for the given set of initial conditions, and cannot 
be continued. However, the inlet mass flux or Mach number can be adjusted 
and iterated if necessary to obtain a flow which is choked at any given 
axial location. In addition, if it is necessary to pass through the throat 
of a choked nozzle, this could be accomplished by local perturbation of the 
throat geometry in a sianner similar to that used by Kreskovsky and Shamroth 
[7] in two dimensions. Finally, if the flow is purely supersonic, 
there is no need for any pressure approximation such as Eq. (2.13), and such 
flows have been treated by McDonald and Briley [8]. 

Arguments favoring subsonic pressure approximation as in Eq. (2.13) lose 
some validity in strongly curved corner regions, where strong cross flows in 
the invlscld outer region of shear layers are deflected. Nevertheless, the 
approximation may not be seriously in error when considering its overall and 
somewhat limited role in the present context. Finally, this t3rpe of pressure 
approximation is questionable for strongly deflected inviscld flows having 
large transverse vortlclty. 

Combining Eqs. (2.10) and (2.13) and setting F ■ F* provides an equation 
nominally governing u^: 

T,*[(U*V)U +(VPy)//> -V(Uj-0j)/2] 

where U is to be written as 

U • •, UjUy + igtVjUy + V,) + I j(WjUy + Wj) 

where u^, v^, are components of U^, and v^, w^ are components of U^. 
Secondary Vortlclty 

An equation governing in compressible flow can be obtained from an 
approximate application of secondary flow theory. The equation governing the 
growth of vortlclty along a streamline for compressible flow with constant 
viscosity is given by Lakshmlnarayana & Horlock [9] and may be written for 
variable viscosity as 


(2.15) 


(2.16) 


8 


(2.17) 


d 


ds 




pq - 


Vx Vp 

P 




where q Is velocity magnitude, s Is the unit vector along a streamline, s Is 
distance along a streamline, 12^ Is vortlcity in the direction of the unit 
principal normal vector n of the streamline, for which R is the principal 
radius of curvature. These quantities are related by the Frenet formula 
n/R = 3s/ds = s • Vs. In Eq. (2.17), G = V x F, and the term containing p 
vanishes if p is constant. Since an intrinsic coordinate system formulation 
as in Eq. (2.17) provides an "inverse" coordinate system if used to compute 
U, and since intrinsic coordinates are degenerate on no-slip surfaces and 
nonorthogonal for general rotational flows, intrinsic coordinates are not 
attractive for numerical computation. However, if the coordinate system used 
for computation is approximately aligned with the flow direction, then an 
approximate equation governing P, can be derived from Eq. (2.17) by replacing 
s by 1, as in the following development; 


Xlg S- iX 

q s ‘U 


ii ' n XI I 

ttd -‘ T 7 


(s.vs).ft (i,-vi,)-n 

qR S -U 


XI 


ni 


U,R, 


(2.18a) 


(2.18b) 


where u, = u_u ; R, is the principal radius of curvature of the x, coordinate, 
II v l ^ 1 

and is vorticity in the direction of n^^, the principal normal of the x^ 
coordinate line. The quantities n^ and R^ are defined by the Frenet formula 
n^/Rj^ *= ij^ • ^i^* illustrate, in an orthogonal coordinate system. 


ax^ h.hj axjJ 


(2.19) 


where h. , h^, h. denote metric coefficients. If p varies, p is replaced by 

^ 2.2 

the imposed pressure p^ as defined in Eq. (2.13). Finally, taking q * 

and s * i^ in the last two terms in Eq. (2,17), and neglecting streamwise 
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diffusion, Eq. (2.17) becomes 


0-7 


( 







« 




+ 



( 2 . 20 ) 


where G* does not contain streamwise diffusion. 

The transverse vortlcity 0^^ In £q. (2.20) contains components which, in 
orthogonal coordinates and assuming an Irrotational are given by 






h, dx, 


I d 

h h- dx, ^*’s'* 8 * 

I d I 


(2.21a) 


A 


3 


jj-yxu 


, Vr 

hg aXg h, dx, 


+ 


I d 

h,hg 7x7^^*''*^ 


(2.21b) 


In the applications contemplated here, the first term on the right->hand side of 
each of Eqs. (2.21a- b) is expected to dominate, and the remaining terms may 
be neglected as a convenience, if desired. 


Energy Equation 

Since the flow being considered is nonadlabatlc, solution of an equation 
governing energy is required. The energy equation can be written in a variety 
of forms, one of which is 

/)0 • VE ■ V- kVT + U • F + <C> ( 2 . 22 ) 

where £ is total enthalpy, T is temperature, k is thermal conductivity, and 
^ is the dissipation function. Solution of Eq. (2.22) by forward marching 
integration requires only that terms representing streamwise conduction of 
heat and also streamwise viscous diffusion in F be neglected. 

Compressibility Relations 

In the imposition of streamwise pressure gradients compressibility 
effects are represented by introducing the perfect gas equation of state 
p “ PRT. No other assumptions are necessary for consideration of compressible 
flow. For moderate subsonic Mach numbers, Invlscid pressure gradients would 
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Ideally be obtained from a compressible potential flow calculation or other- 
wise from an Incompressible potential flow corrected for compressibility 
using either the Prandtl-Glauert formula or Laltone's variant [10], 

Replacing p in Eq. (2.13) by Che state equation and eliminating temperature 
using the temperature-enthalpy relation 


E 


CpT + 


UU 

T" 


(2-23) 


where c^ denotes specific heat, the following auxiliary equation relating 
the Imposed pressure gradients with density, velocity, and total enthalpy is 
obtained: 

V'^[pl + Pv<x,j] ■ 

where y is specific heat ratio. A slight simplification results if U • U is 
- — 2 

replaced by (i^ * U) in Eq. (2.24). 

Typically, the treatment of nonconstant total enthalpy is of Interest 
for predicting heat transfer levels in shear layers. However, it should be 
noted that the geometrical pattern of streamlines for steady inviscid flow 
of a perfect gas depends only on the total pressure gradient and not on the 
distribution of total enthalpy (cf. Hawthorne [6]). Thus, within an assximp- 
tion of small transverse total pressure gradient in the free stream or core 
region, it is possible to consider core flows with two separate streams having 
large differences in total enthalpy and velocity and hence separated by 
thermal and viscous shear layers, provided the two streams have about the 
same static and total pressure. In addition, the analysis can be applied 
to flows with two streams at different total pressure provided the static 
pressure is continuous across the Interface. Although Inviscid pressure 
gradients can no longer be represented by a two-dimensional potential flow, 
the error thus Introduced by the present analysis may not be serious. Since 
the mean pressure level is determined accurately by p^, the primary source 
of error in this Instance would be that due to the assumed transverse variation 
in imposed streamwlse pressure gradients, which in turn is associated with 
streamline curvature. For geometries without appreciable turning, this error 
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■ay be acceptable. The foregoing observations regarding two stream flows were 
used to advantage in the present application to flow in turbofan lobe mixers. 

Viscous and Heat Conduction Terns 

A variety of techniques for suppressing streaowise diffusion is possible, 
and- the particular choice in any given application may depend on the complexity 
one- is willing to tolerate in the viscous terms. For example, the complete 
set of viscous terms can be written out and all derivatives in the streamwlse 
direction neglected. Alternatively, only second derivatives in the streamwlse 
direction may be neglected. Here, attention is restricted to incompressible 
flow but with variable viscosity. In the momentum equation (2.10), 

F s -V X vV X U , where v is kinematic viscosity. A particularly simple 
expression which neglects contributions from is given by 

t,' F * - I,* VXvVxtjUjUy • P' (2.25) 

Alternatively, contributions from may be retained as in the approximation 

f,»P * vVjXU ■ P' (2.26) 

Similarly, in the vorticity equation (2.17), 6 is defined by G « -V x V x vO 
and approximated by 

Ij.g i -f, *7X7X1, -G (2.27) 

Streamwlse heat conduction in Eq. (2.22) is suppressed by the approximation 

7- k7T * 7*k7,T (2.28) 


Governing System of Equations 

A complete system of six coupled equations governing u^, (2^, «, E, 
and p is given by Eqs. (2.8), (2.9), (2.15), (2.22), (2.24), and (2.20). 
Ancillary relations are given by Eq. (2.5) for composite velocity, Eq. (2.14) 
for mass flux, and Eqs. (2.21a-b) for transverse vorticity. 
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The Turbulence Model 


The present analysis is to be used to predict a turbulent flow field, 
thus the stress terms F* and G' In Eqs. (2.15), (2.20), and (2.22) contain 
turbulent shear terms in the form of fluctuating velocity and velocity- 
temperature correlations. These terms may be modeled using any of a number 
of turbulence models. In the present application to lobe mixer flows, these 
turbulent shear terms are modeled through the introduction of a turbulent 
eddy viscosity determined from solution of two partial differential equations 
governing the transport of the turbulence kinetic energy and the turbulent 
dissipation. The particular model used in the present application is described 
by Launder and Spalding [11]. Using Cartesian tensor notation, and overbars 
to indicate time averaging the turbulence kinetic energy is defined as 


and the dissipation as 


k = 


2 

2 


Uj Uj 


€=V 


du\ du'j 
dx j dxj 


(2.29) 


(2.30) 


By hypothesis, Prandtl and Kolmogorov have suggested that the effective 
viscosity is proportional to the local density and the product of a character- 
istic turbulent velocity and length scale. The characteristic velocity is 

assumed equal to the square root of the local value of turbulence kinetic 

1/2 

energy, thus = cpJik where c is a constant of proportionality. Through 
dimensional arguments, the length scale may be related to the turbulence kinetic 
energy and dissipation as 


C 




(2.31) 


and from the Prandtl-Kolmogorov hypothesis it follows that the turbulent 
viscosity may be expressed in terms of k and c as 


C*pk^/t 


(2.32) 


Based on an examination of a large amount of experimental data Launder [11] has 

4 

recommended that c ^ 0.09 = c , and under these conditions the length scale 

V 

may be thought of as the conventional mixing length. 
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The transport equations governing k and e at high Reynolds number are given 
in vector form £L1] as 


P—’¥ pU *Vk* Vk +P- pc (2.33) 

d\ 

p ^ +fU'7«i7-(/i+/i^/C’j)VK + C * P-Cj^— (2.34) 

il k k 


In Eqs. (2.33, 2.34), P is the turbulence production and is given in the APPENDIX. 

The quantities , Ci, C 2 are empirical constants which, as recommended by Launder 

& Spalding [ll] , are taken as 1.0, 1.3, 1.44, and 1.92, respectively. 

Although the k-£ turbulence model provides a general means by which turbulence 
effects may be modeled, experience with the model has indicated that it may not be 
completely reliable. Negative values of dissipation may develop under certain con- 
ditions at isolated points in the flow field leading to a breakdown of the calcula- 
tion. To provide a means to circumvent this problem an optional turbulence model 
was included in the analysis. 

This turbulence model is referred to as a wake turbulence model and is derived 
from the definition of turbulent viscosity Pt Prandtl's suggestion that the 
dif fusional flux of some property ^ ^ ' may be expressed 



I ,1/2, ^ a<^ 

By ' P ~ST 


hence one obtains the Prandtl-Kolmogorov relationship 


' k 

9 


(2.35) 


(2.36) 


The wake turbulence model is obtained by assuming that the turbulent length scale 
and velocity scales are constant throughout the flow field. From equations 2.31 
and 2.32 it follows that the constant 1/a^ is simply The velocity scale is 

taken as the square root of the initially specified turbulence kinetic energy k. 
The length scale is specified using the geometric constraints of the mixer as a 
guide, for instance assuming that the initial length scale would not be smaller 
than the inlet boundary layer or larger than the duct height. As a consequence 
of this formulation, the turbulent viscosity is dependent only upon the initial 
choices of k and I , and the local density which is determined as the solution 
progresses. 
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APPLICATION TO FLOW IN A LOBE MIXER 


In this section, specific details arising In the application of the 
foregoing analysis to lobe mixer problems are given. A cross section of a 
typical lobe mixer geometry Is shown In Fig. 2. The area Immediately 
downstream of the nozzle plug Is faired In with an assumed streamline to 
model the reversed flow region expected In this region. Since the flow area 
thus excluded from consideration is very small, this treatment Is not believed 
to ^xitroduce significant error. A more detailed analysis of this separated 
flow region could be performed following a zone embedding approach using the 
Navler-Stokes equations (cf. [12]). Although the mixer geometry Is axlsymme- 
trlc, the flow Is three-dimensional due to the azimuthal variation of the hot 
and cold streams. However, due to observed symmetry, only a pie-shaped region 
of the transverse coordinate surface need be considered. The shape of this 
region and the extent of typical hot and cold streams at the mixer exit 
surface are shown in Fig. 3. 

Coordinate System 

Curvilinear orthogonal coordinates x, y, z are constructed to fit the 
flow passage boundaries as shown in Fig. 2. Metric coefficients h]^, h 2 » h 3 

are defined such that Incremental distance s is determined by (6x)2 « (h. 6x)^ 

2 2 ^ 

+ (h^^y) -f (h^5z) . In planes of constant z, azimuth, orthogonal streamlines, 

and velocity potential lines from a two-dimensional planar Incompressible 

potential flow analysis are utilized as the coordinate lines for constant 

y and x, respectively. This x-y coordinate system is then rotated about an 

axis representing the centerline of the lobe mixer, to form an axls 3 nmnetric 

coordinate system. Thus the z direction can be regarded as cylindrical 

(l.e., 6) and h^ ■= r. The x coordinate Is taken as the primary flow or axial 

coordinate, and Is associated with surfaces for which the two-dimensional 

potential is constant. Ttie y and z coordinates define transverse secondary 

flow planes at any given x location. In this coordinate system, the normals 

1 to the transverse planes coincide with 1. the unit vector In the x direction, 
n j. 

A two-dimensional incompressible potential flow analysis and computer program 
developed by Anderson [13] was employed In the present Investigation, without 
modification, to compute the necessary coordinate data. 


15 


The Axlsynmetrlc Pressure Field 


Since the coordlnete system is based on a two-dlnenslonal planar 
potential calculation, the coordinate lines are not precisely aligned with 
the axisyBBietrlc potential flow appropriate for the alxer geometry. For the 
purpose of setting the primary and secondary flow directions, however, the 
two-dimensional solution provides a good approximation and thus the direction 
of ip the unit vector in the x coordinate direction, is taken as the primary 
flow direction. This assumption is Justified since surfaces of constant 
potential from the two-dimensional solution remain surfaces of nearly 
constant potential even though the value of the potential at each surface may 
be altered significantly in the axlsymmetric case. The two-dimensional 
potential flow is not suitable for approximating the axlsymmetric potential 
flow pressure gradients, however, and consequently an axlsymmetric potential 
flow is confuted in the given coordinate system, prior to the forward 
marching calculation. 

Since a compressible potential flow solver was not readily available, 
an incompressible axlsymmetric potential flow was computed and then corrected 
for compressibility during the forward-marching solution procedure, using 
laitone's rule (Ref. 10). Although the Laitone compressibility correction is 
expected to become Inaccurate at high subsonic Mach numbers, this error is 
mitigated both by the computed correction p^ for mean pressure drop, and by 
certain mixer flow characteristics. Since the mean pressure level is deter- 
mined accurately by p^ at each axial location, influence of the Imposed 
inviscid pressure p^ is limited to approximating the transverse variation of 
axial pressure gradient and hence is associated only with curvature of 
potential flow streamlines. The axial development of mean pressure drop and 
average ^lach number is thus Influenced only indirectly by the Imposed pressure 
and should be predicted with reasonably accuracy. Furthermore, in siost mixer 
geometries, curvature effects are more pronounced in the upstream region and 
near the end of the center plug, i^ere the Mach number is relatively low. 

In the region of strong acceleration and high subsonic Mach number near the 
exit, streamlines tend to have small curvature, and hence the flow has small 
transverse variation in axial pressure gradients. 

The incompressible potential flow is determined by solution of the 
axisymnetrlc potential equation 
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(3.1) 


using an efficient Iterative ADI procedure and the solution, which represents 
the correct axlsynimetrlc potential, Is used to obtain the axlsymmetrlc pres~ 
sure field, P^(x,y). In setting boundary conditions for the axisyunnetric 
potential, the Inlet and exit planes are assumed to be surfaces of constant 
(but differing) potential, and normal derivatives are set to zero elsewhere. 

Since £q. (3.1) may be considered a dimensionless equation with the 
dimensionless potential given by 


<k 


<#-* 


<3. 2) 


* 2 / 

where <p has units of ft /sec, the resulting potential field may be arbitrarily 
scaled to give a mean dimensionless velocity of unity at the mixer exit plane. 
This normalized potential field is then used in the computation of P^(x,y). 


Governing Equations in Orthogonal Coordinates 

A general model for three-dimensional flows which may have large secondary 
velocities was given in vector form in the preceding section. The governing 
equations are given here for the orthogonal coordinate system used to 
represent the mixer passage geometry. Throughout the remaining discussion, 
all variables in the governing equations are nondimensional, having been 
normalized by the following reference quantities: distance, Lj.; velocity, Uj.; 
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density, p^; temperature, T^; total enthalpy, U^; pressure, viscosity. 

Here the subscript r denotes a reference quantity. This normalization 
leads to the following nondimensional parameters: Mach number, M; Reynolds 

number. Re; Prandtl number, Pr; and specific heat ratio, y. These parameters 
are defined by 

M . U,/C , Re , Pr - Cp/i,/k ,yCp/ (3.3) 


where is the molecular viscosity, k is thermal conductivity, and c^ and c^ 

are the specific heats at constant pressure and volume. The reference speed 

2 

of sound, c, is defined by c yR T , where R is the gas constant. Since 

R t g 

the flow is assumed to be turbulent, the dependent variables are taken to be 
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the tlBie— everaged quantities in the usual sense j; however, third order correla- 
tion as well as second order correlations Including the fluctuating density 
have been neglected. The decomposition of the velocity vector thus can be 
written 

0 ■ IjU + 7jV + TjW (3.4) 

Where T^, and are the unit vectors in the x, y, and z coordinate 
directions, respectively. Since is closely aligned with the primary flow 
direction, v^ is small and is neglected as a minor convenience. This 
approximation is of little consequence since the neglected quantity v^ 
effectively reappears as part of v^ and is thus determined from solution of 
the continuity equation during the forward— marching solution process. Instead 
of being imposed _a priori from the potential flow. Since Wj is zero by 
definition, the velocity decomposition can be written as 

U -TiUp+itV, +T,w, (3-5) 

Since u is the only streamwise component of U the subscript "p" is 
P 

dropped. Under the stated assumptions the streamwise momentum equation can 
be written as 


du du du 


dh, 2 dh 
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II 


The energy equation is approximated by 
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(3.7) 


with the Prandtl number Pr and turbulent Prandtl number Pr^ assumed to be 1.0. 
From the axisymmetrlc potential flow. 


I a<^ \ dcf> 

“i “ “h^ *ajr * ''i * "h^ 


(3.8) 


Substituting these velocity components into the inviscid, incompressible 
momentum equation 


/ Ut-U, \ _ - . • „ 

v(-s^)-u^x{vxup-- — VP„^, 


(3.9) 


and integrating Eq. (3.9) with the observation that V x Uj = 0 the inviscid 
incompressible pressure field Pj^j_p(x,y) is obtained from the relationship 


^(u/+Vj2) 

^ ' constant - 


(3.10) 


The imposed compressible pressure gradients are then determined from P.. ^ as 

M*U 

i * ^^M»0 (3.11) 

dx yi-M2 dx 

where is the local Mach number. The Imposed pressure gradient is related 
to the dependent variables using the following form of the auxiliary gas law 
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dP. 


I ^ 


ax 


dx 


(3.12) 


The integral Bass flux relation is 

ff hjhs^udydz »C 

To obtain the rotational secondary flow components of the velocity, the 
streamwise vortlclty equation is written as 

U/ d€ ^dpu\ . V, h, I af ^ dpu \ 

0 ( dx ^ dx / o hj V dy dy / 

”o” hj { az ^ az / * hj ay 


(3.13) 
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where Q - u **" vorticity is related to velocity by 
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The vector potential ij/ is determined 

from 


1 « 1/ »»» \ 1 
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ay J 



(3.15) 


-CC3-16) 


and the rotational secondary flow velocities are determined as 

" 7*^ ~dT“ 


(3.17a) 
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(3.17b) 


The scalar potential is governed by 
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(3.18) 

dx 


(3.19a) 


(3.19b) 


The secondary flow velocities are thus 

V* - 

*is ‘ ”<#> 


(3.20a) 

(3.20b) 


The turbulence model equations (Eqs. 2.33 & 2.34) must also be expressed in 
orthogonal coordinates. It is assumed that the turbulence equations are 
tensor invariant [14] and thus may be expressed directly in orthogonal 
coordinates by performing vector operations in such a coordinate system. 
With streamwise diffusion neglected to allow forward-marching solution, the 
result for the turbulence kinetic energy equation in steady flow is 


+ 


dk dk dk 



(3.21) 
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and almilarly the turbulence dleslpation equation is given by 




1 

— 1 




Re 

dy 1 
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(3.22) 
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- hjhjhjCgP ^ + C| ^ p 

The turbulent viscosity is then 

C^^k* 

- -^R. 


(3.23) 


and the length scale may be determined as 

3M.3/2 
C, 

I 


Cfj, k 


(3.24) 


Numerical Method 

The governing equations are approximated by finite differences. Three 
point central differences are used for all transverse derivatives. An 
analytical coordinate transformation devised by Roberts [15] can be employed 
as a means of introducing a nonuniform grid in each transverse coordinate 
direction to concentrate grid points near vails while maintaining second- 
order accuracy. In the computation presented here a uniform mesh spacing 
in each transverse direction was used. The streamvlse vorticlty equation 
is decoupled from other equations in the system and linearized with respect 
to ( by an ad hoc process consisting of lagging quantities not yet available 
at the implicit level. The resulting implicit difference equation for C is 
solved using a scalar ADI scheme based on the technique of Douglas & Gunn [16] 
for generating ADI schemes as perturbations of fundamental implicit schemes. 
Equations (3.16 & 3.18) for the vector and scalar potential functions ^ and 
are elliptic in the transverse planes and are solved, given values for 
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the right-hand sides, using scalar Iterative ADI. Specifically, the Douglas- 
Gunn perturbation of the Crank-Nicolson scheme is used for Eqs. (3.13 & 3.16) 
while perturbation of the backward difference scheme is used for Eq. (3.18). 
With regard to the streamwise momentum and energy equations (3.6 '& 3.7), since 
the transverse velocity components are determined separately from elliptic 
equations the coefficients of the transverse convective terms are lagged. 

Formal linearization of the streamwise convective terms following Briley and 
McDonald [17] as extended to nonconservation law form by Kreskovsky and 
Shamroth [7] also indicates that to first order the coefficients of these 
derivatives may also be lagged. The energy and momentimi equations thus 
become decoupled while maintaining first-order accuracy in the spatial marching 
direction. The resulting difference equations may then be grouped by coordinate 
direction and solved using an efficient ADI scheme. The turbulence model 
equations form a similar system, but these equations remain coupled through 
the source terms. In this case, the resulting difference equations can be 
written in block-tridiagonal matrix form and solved efficiently using linear- 
ized block implicit (LBI) techniques. 

A summary of the overall algorithm used to advance the solution a single 
axial step follows. It is assumed that the solution is known at the n level 
x^ and is desired at 

1) The imposed pressure gradients are determined from Eqs. (3.10) and (3.11). 

2) A value for the mean viscous pressure drop dp^(x)/dx is assumed. 

Initially, the value from the previous step is used. 

3) Equations (3.6), (3.7) and (3.12) are solved to obtain values of u^"*"^, 

E^"*"^ and In general, the integral mass flux relation (3.13) will not 

be satisfied. 

4) Return to step (2) and repeat this process iteratively using the standard 

secant method [18] to find the value of the mean viscous pressure drop which 
leads to u*^^^ and which satisfies the integral mass flux relation (3.13). 

The secant method was found in practice to converge the integral mass flux to 
five figures on the third iteration, unless the flow is near transonic. 

5) The turbulence model equations (3.21 & 3.22) are solved using a split 
LBI scheme. 

6) The streamwise vorticity equation (3.14) is solved using scalar ADI to 
obtain 
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7) The vector potentiel equation (3.16) is solved using scalar iterative 
ADI to obtain 

8) Values for solenoldal secondary flow velocities v^ and are conputed 
using Eq. (3.17). 

9) Using values now available for 8pu/3x» the scalar potential (3.18) is 

solved using scalar ADI to obtain to assure that the solution at level 

n4>l satisfies continuity on a differential basis, and the irrotational 

secondary flow velocity components v. and w. are determined from Eq. (3.19). 

9 9 

To ensure that the initial conditions satisfy continuity on a differential 
basis, steps 2 through A and step 9 are performed to obtain at the initial 
station either using estimated values of the transverse velocities (obtained 
from experimental data) or assuming the transverse velocity components are 
aero. 

Specification of Initial Conditions 

The initial conditions for a lobe siixer calculation may be specified by 
using either an automated starting routine or by reading in experimental 
data. The fully automated procedure is considered first. 

To obtain initial conditions for a lobe mixer calculation, it is necessary 
to specify the velocity and total temperature of the respective hot and cold 
streams in addition to a mean value of static pressure. The lobe shape is 
specified, and based on the specified lobe shape, a decision is made as to 
idiether a grid point lies in the hot or cold stream. The appropriate values 
for velocity, temperature and density are then assigned to the grid point. 

The total energy is assumed to be constant but may differ in both streams. 

The velocity in the two respective streams is assigned its nondimenslonal 
reference value which is then corrected to account for normal pressure 
gradients present at the initial plane as determined from the axisymmetric 
potential flow. To account for boundary layers on the lobe, hub and shroud 
surfaces, the free stream velocity profiles are further scaled in accordance 
with an assumed turbulent boundary layer velocity profile, and the distance 
from the lobe surface. 

Referring to Fig. 3, the position of the line representing the lobe 
shape, S- , is determined by specifying its r and 6 coordinates at a sufficient 
number of points to describe the lobe shape, thus 
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- s^(r,0) 


(3.25) 


. . or along 9 
min ® max 


max 


where S may have an arbitrary shape subject only to the following restrictions 

JLi 

(see Fig. 3) . 

1) must enter the computation region along r^ 

2) Sj^ must exit the computation region along 6 ^^ ^ or along r^ 

3) In the specification of the r, 0 coordinates of S , r must increase 
monotonically, i.e., any line of constant radius may intersect S_ only once 

Li 

when moving along S in the direction indicated in Fig. 3. No restriction is 

Li 

made on e. An example of the lobe shapes which may be specified is shown in 
Fig. A. The initial velocity profile with boundary layers are given by 


u(u 


STRM 


• V-y/s., ,s,H) 


(3.26) 


where is the velocity of the appropriate stream, is the boundary 

oLKM fiL 

layer thickness, y/6 represents the shortest distance from a surface (hub, 

it 

shroud or lobe) to the grid point in question, 5 is the local displacement 
thickness and H is the shape factor. For grid points near the lobe the 
distance, y/5, is determined as shown in Fig. 3a by constructing a triangle 
given two data points on the lobe and the grid point in question. The normal 
to the base of the triangle (the line segment between the two lobe data points) 
is constructed and its length determined. If a normal cannot be constructed 
such that it Intersects the line of which the base is a segment between the 
lobe data points, the distance to the lobe is taken as the distance to the 
nearest data point as shown in Fig. 3b. For grid points near the hub or shroud 
the calculation of y/S is straightforward. Thus having determined y/6, the 
velocity at the grid point in question can be scaled to account for the 
boundary layer. The assumed boundary layer profile used in the present code 
is a Coles-type profile modified as suggested by Waltz [19]. 

u - u^[(-.3y*“ 5)e”®^ ^ (3.27) 


Since the turbulence model used is valid only at high Reynolds numbers, the 
low Reynolds number sublayer flow in the immediate region of walls must be 
treated In an approximate manner. This is accomplished by using wall functions 
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which assume that locally between the wall and the first grid point away 
from the wall the velocity profile is logarithmic such that 

I . 




(3.28) 


where is the friction velocity, k ■ 0.43 is the von Karman constant, and 
c is taken as 5.0. It follows that 


du 

ly 


Ky 


(3.29) 
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(3.30) 


Eqiiation (3.29) is used to determine a wall slip velocity such that the finite 
differenced form of the velocity gradient one point off the wall is consistent 
with the assumed law of the wall profiles. The density is then computed to 
be consistent with the temperature and the local static pressure as determined 
from the mean static pressure corrected for the normal pressure gradients. 

Initial conditions for the turbulence model are determined through speci- 
fication of a reference length scale, And free stream turbulence levels. 

Near the hub and shroud the length scale is determined from the McDonald- 
Camarata [20] relationship 

IIL 

I 

'ref 

where y is the distance from the surface in question. In the region about the 
lobe surface the distribution of Eq. (3.30) may be used with y representing 
the distance to the lobe from the grid point in question. However, since the 
lobe may pass arbitrarily close to a grid point the length scale determined 
from Eq. (3.30) will approach 0.0 as y approaches 0.0. This results in a low 
Reynolds number region in the initial flow field which can create problems 
with the turbulence model. This problem is circumvented through use of an 
option which renders the length scale constant in the region near the lobe, 
and the distribution of Eq. (3.30) is used only near the hub or casing if 
boundary layers are specified there. With the length scale distribution known, 
the initial turbulent viscosity is obtained from the generalized mixing length 
relationship 

1/2 

<■,■^<■[(•5; 0 * * (if?) 


26 



The initial turbulence kinetic energy and dissipation are then determined from 
simultaneous solution of Eqs. (3.23) and (3.24). The kinetic energy is 
modified slightly following McDonald and Kreskovsky [21] to account for the 
free stream turbulence, and the turbulent viscosity is then recomputed from 
k and e to be consistent. 

Under certain conditions, neither the wall length scale option £q. (3.30) 
nor the wake length scale option (constant length scale) will provide satis- 
factory initial values of k, e and ji^. This may occur when the mesh is coarse 
and the velocity gradients in Eq. (3.31) vary greatly between adjacent grid 
points. Under these circumstances, it can be useful to initiate the turbulence 
quantities using an option which specifies the turbulent kinetic energy, 
dissipation and viscosity as essentially constant throughout the flow field. 

Such an approach has also been used previously by Launder, Morse, Rod! and 
Spalding [22] in the computation of free shear flows. 

If experimental profiles are read in, the specification of initial 
conditions is similar. The velocity and static temperature are input at each 
grid point. The lobe shape is specified as in the automated starting proce- 
dure, and a reference length scale and free stream turbulence intensities must 

be provided. Additionally, estimates of the boundary layer thickness on the 

% 

hub, shroud and lobe must be specified for use in setting up the initial 
turbulence quantities. Initialization of the turbulence quantities is similar 
to that used in the automated procedure. 

Boundary Conditions 

To march the solution downstream, it is necessary to apply boundary 

conditions in each transverse plane. Boundary conditions are needed for the 

energy equation, the primary flow momentum equation, the streamwise vorticity 

equation, the vector and scalar potential equations, and the turbulence model 

equations. Referring to Fig. 3, the boundaries at 6 . and 6 represent 

mxn max 

symmetry planes. Boundary conditions on these surfaces are therefore straight- 
forward and are as follows: 

1) azimuthal gradients of u , E, k and e are set to zero, 

2) the vorticity is set to zero, and 

3) the normal velocity component is zero, and thus the vector potential is zero 
and the azimuthal gradient of the scalar potential is zero. 
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Boundary conditions on the hub and shroud surface are specified as follows: 

1) the nortnal velocity is sero along these surfaces, and the vector potential 
and the normal gradient of the scalar potential are thus set to zero, and 

2) these surfaces are taken as adiabatic, and thus the normal gradient of 
energy is set to zero. (Although not included in the present version of the 
code, thermal vail functions can be developed which would allow specification 
of the wall temperature. These boundary conditions could be easily 
implemented.) 

Boundary conditions for the primary flow velocity, u , and the turbulence 
quantities k and e are treated in a manner consistent with the high Reynolds 
number limitation of the turbulence model. These boundary conditions are 
based on the same logarithmic profile assumption used in the Initial conditions 
(Eq. (3.28)). The gradient condition of Eq. (3.29) is used to determine a slip 
value of the primary velocity, u^. Furthermore, consistency with the log law 
Implies that the turbulence is in local equilibrium and that the shear stress 
is constant near the wall. From this it may be deduced that 


in the near wall region, and 


C 1/2 
> 


(3.32) 


dk _ (3.33) 

-r — ■ O 

dy 


Equation (3.32) is used to specify the value of k one point off a wall with 

Eq. (3.33) giving the wall value. Additionally, the log law implies that the 

length scale, £, varies linearly with distance from a wall (£ * <y) . Thus 

initially the turbulent dissipation one point off a wall was specified using 

Eq. (3.24). However, problems were encountered with this boundary condition 

which led to negative values of e. This problem was traced to the linear Iza- 
3/2 

tion of k in Eq. (3.24) and was resolved by combining Eqs. (3.24) and (3.32) 
to give the linear relationship 

U ’ (3.34) 

* “ Ky 

which performed satisfactorily. 

With regard to the boundary conditions for the secondary vorticity, it is 
observed that the generation of vorticity is most likely to occur in the free 
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shear layer downstream of the lobe rather than within the hub and shroud 
boundary layer. Thus the secondary vorticity is assumed zero in the near wall 
region. It should also be noted that as the flow proceeds past the hub 
region, provisions are made to change from wall boundary conditions on the 
hub to symmetry conditions at the centerline. 
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COMPUTED RESULTS 


Several flow calculations were perfomed to assess the computational 
procedure and turbulence XK>del. These flows Include a simple coaxial 
turbulent Jet, used to validate the turbulence model, and both laminar and 
turbulent mixer calculations In two different geometries typical of modem 
designs. 

Coaxial Jet Flow 

The first flow considered Is that of a coaxial Jet experimentally 
Investigated by Forstall and Shapiro [23] and used as a test case In the NASA 
Langley Free Turbulent Shear Flow Conference [24]. This test case was 
chosen In part since other predictors at this conference were able to obtain 
good agreement with this data, using the same two-equation turbulence model 
as Is used in the present lobe mixer program. These experimental measurements 
were obtained with a .25 In. diameter nozzle exhausting into a 4 In. diameter 
pipe. The velocity ratio of the nozzle stream to the outer stream was 4 to 1. 
Tabulated values for the initial streamwlse velocity profile are given in 
[24] as are values for the Jet spread rate and the centerline velocity decay. 
Profile data at various dowstream locations may be found In [23]. 

Since the shear layer is very thin in the upstream region near the Jet, 
a very fine mesh would be required to resolve the flow both near the jet and 
at large distances downstream, where the shear layer thickness Increases 
many fold. For convenience, the flow was computed in two parts with different 
outer boundaries and grids, upstream and downstream of X/D • 1.0. The first 
calculation had an outer pipe radius of 16 nozzle radii, utilized 20 radial 
grid points concentrated near zr/R^ “ 1, and was started at X/D ■> 1. The initial 
velocity profile was chosen to agree with the measured profile of Sami, Carmody 
and Rouse [25] at X/D ■ 1.0. Turbulence quantities were initialized using a 
constant length scale across the shear layer In conjunction with Eqs. (3.23), 
(3.24) and (3.31). The calculation was marched to X/D “ 37. 

Comparisons of the predicted and measured centerline velocity decay and 
jet spread rate are presented in Figs. 6 and 7, respectively. The predicted 
centerline velocity decay is in very good agreement with the measurements, 
and the predicted spread rate Is only slightly above the trend of the data. 

It should be noted that these calculations did not utilize the modifications 
to C 2 and c^ suggested by Launder and Spalding [9] to Improve predictions of 
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the spread rate. Recently, Janicka and Kollmann [26] also observed good 
agreement between predicted and measured spread rate and velocity decay using 
the unmodified turbulence model. Figure 8 shows the predicted streamwise 
velocity profiles Cnormalized by jet half radius) at various X/D locations 
and the range of data measured by Forstall. Again the agreement is good 
despite a slight under resolution of the outer region at X/D * 37, which is 
evident from the grid point locations in the figure. 

The second part of this flow calculation was performed in the upstream 
region 0 i X/D i 1.0, with a computational domain limited to 0.5 i 
As shown in Fig. 9, the initial profile was chosen to match the experimental 
data of Forstall and Shapiro [23 & 24]., and includes a small slip velocity 
at the nozzle edge. The distribution of grid points is shown in Fig. 9 
and provides adequate resolution in the region near the nozzle. Turbulence 
variables were again initialized using a constant length scale. In Fig. 10, 
the computed velocity profile at X/D = 1.0 is compared with that used as the 
initial profile in the previous calculation with expanded computational 
domain. The good agreement between the velocity profiles in Fig. 10 justifies 
the initial conditions used for the previous solution obtained in the region 
downstream of X/D = 1.0. The calculations of coaxial jet flow given in this 
section confirm that the turbulence model being used perfoims satisfactorily 
for this relatively simple flow for which adequate experimental documentation 
is available. Further, the predictions of the lobe mixer deck for this case 
are in essential agreement with predictions made by other investigators using 
the same turbulence model [24]. 

Lobe Mixer Flows 

Several calculations of lobe mixer flows were performed to demonstrate 
the overall method in its present state of development and to explore the 
potential of the method for making detailed predictions of the flow and mixing 
process. A stimmary of the calculations performed is presented in Table I. 

The test cases include both hot and cold flow calculations in two different 
geometries, and with different inlet velocity ratios and average Mach number. 
All calculations were made using 20 equally spaced grid points in the radial 
coordinate direction and 10 equally spaced points in the azimuthal direction. 
Twenty to twenty-two axial stations were computed. 
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Before discussing results fron these calculations. It is useful to 
define average values of total pressure loss coefficient, ideal thrust 
coefficient, and Mach number at each axial location, to aid in the evaluation 
of flow properties. A mass-averaged total pressure loss coefficient is 
defined as 




Vef 


pU • dA 


'PT 


J />U*dA 


(4.1) 


where is a reference total pressure associated with the turbine stream 

■^ref 

reference conditions, and snd the reference density and velocity 

associated with the turbine stream. The value of at the starting plane 
is usually not zero due to the presence of shear layers in the initial condi- 
tions, and for clarity this starting value is subtracted from Cp^ in the 
presentation of results. The distributions of Cp^ presented thus represent 
losses Incurred downstream of the starting plane. An ideal thrust coefficient 
is defined and based on the thrust which would be obtained by isentropic 
expansion from local conditions to a predetermined exit pressure. The ratio 
of the local value of this thrust T to the value T^ based on the initial flow 
field is defined as the ideal thrust coefficient, T/T^. The thrust is 
obtained from the relationship 

T - / yP,M,*<JA, (4.2) 

where is the assumed exit pressure, is the Mach number based on 
isentropic expansion from the local pressure to F , and A is the effective 
exit area, again based on isentropic expansion to P^. Finally, an area- 
averaged Mach number is defined as 



JRe^u^t^ 

The first two solutions were obtained for geometry A, which is shown in 
Fig. 11. The lobe shape specified for this geometry represents a high bypass 
ratio configuration and is shown in Fig. 12. As indicated in Table X, Case 1 
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is a laminar flow calculation at a Reynolds number of 1280 » based on fan 
stream properties and the Initial mixer diameter. Boundary layers on the hub 
and shroud were neglected In this calculation, and Instead slip conditions 
were imposed on these boundaries. Case 2 Is identical to Case 1 except that 
the flow is turbulent flow with Reynolds number of 2.9 x 10^, and plug, 
casing and lobe boundary layers are Included. In each case, the ratio of 
fan to turbine temperature is 0.5, and the ratio of fan to turbine velocity 
is 0.707. These conditions result in matched Mach numbers for the fan and 
turbine streams. The average inlet Mach number is 0.34 in each instance. 

The streamwlse distribution of total pressure loss coefficient is shown 
in Fig. 13 for Cases 1 and 2. The laminar flow has higher loss levels despite 
its neglect of wall boundary layers, a consequence of the lower Reynolds 
number. The predicted axial variation of the ideal thrust coefficient is 
shown in Fig. 14. Here it is observed that the turbulent flow calculation 
exhibits a significantly higher gain in gross thrust than the laminar calcula- 
tion. These higher thrust levels prestunably reflect both increased mixing 
and lower losses in the turbulent flow case. It is also noted that the 
exit Mach number for these two cases was about 0.5 to 0.6, whereas actual 
full scale mixers would normally be choked at design conditions. 

Re^u^t^ _f^o^ 

The remaining mixer flow calculations were performed for geometry B, 
which is shown in Fig. 15. Two lobe shapes were considered, and these are 
shown in Fig. 16. Each of these lobe shapes provides a mixer with high lobe 
penetration and a bypass ratio significantly lower than that of geometry A. 

The first set of flow calculations considered here consists of Cases 3 
and 4, which correspond to lobe shape B-1 in Fig. 16. These two cases 
represent both hot and cold flow conditions for which the inlet Mach numbers 
(and total pressures) of the fan and turbine streams are matched. In each 
case, the average initial Mach number is 0.194 and was estimated to give 
nearly choked flow at the nozzle exit. The fan stream Reynolds number is 
2.6x10^ in each case, and boundary layers on the plug, casing and lobe were 
included in all of the turbulent flow calculations. 'Test Case 3 is typical 
of flight cruise conditions and has a fan to turbine temperature ratio of 0.4 
and velocity ratio of 0.7. In Case 4, the fan stream conditions were identical 
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to those of Case 3, but the turbine stream temperature was reduced to match 
that of the fan stream, and in addition the velocity of the turbine stream 
vas adjusted to produce the same Mach number and total pressure ratio as 
Case 3. Taken as a set. Cases 3 and 4 simulate a common experimental procedure 
for testing lobe mixer designs whereby both hot and cold flow measurements 
are taken with the same conditions for stagnation pressure. Since the cold 
flow measurements provide an indication of losses associated with plug, 
casing, and lobe boundary layers, the differences between hot and cold flow 
measurements provide an indication of thrust gains and total pressure losses 
due to the mixing process Itself. The set of test Cases 3 and 4 thus demon- 
strate one potential use of the present calculation procedure for aiding in 
the design and evaluation of lobe mixers . 

Comparisons of the computed axial development of total pressure losses, 

Mach number, and thrust are given in Figs. 17-19, respectively, for Cases 3 

and 4. The total pressure losses for the hot flow case are somewhat higher 

than those of the cold flow case, although the distributions are similar. 

The predicted loss levels within the mixing duct are of the order of one 

dynamic head at fan stream inlet conditions. The higher loss levels present 

in the hot flow case are presumably a result of higher losses in the mixing 

layer between the hot and cold streams, since the difference in velocity 

between the two streams is much higher in the hot flow case. The axial 

variation in Mach number is shown in Fig. 18. Although the initial development 

of Mach number is very nearly the same for the two cases, as a result of 

higher losses the hot flow case undergoes stronger acceleration near the 

exit and attains an exit Mach number of 0.91, as compared with 0.8 for the 

cold flow case. The rapid acceleration downstream of X/L^ « 1.6 occurs 

primarily because of the change in area of the geometry but is due in part to 

frictional losses. On the other hand, since loss levels Increase with Mach 

number, the total pressure losses in Fig. 17 are seen to Increase rapidly 

downstream of X/L " 1.6. 

r 

The final flow parameter of Interest here is the thrust coefficient, 
which is shown in Fig. 19. As is to be expected, the thrust coefficient for 
the cold flow case undergoes a monotonic decrease, as viscous losses accumulate. 
In the hot flow case, the effects of thermal mixing provide a significant 
increase in thrust over most of the mixing duct. Beginning at about X/L^ ■ 2.0, 
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the high loss levels associated with the nozzle-induced acceleration outweigh 
any further thrust augmentation due to mixing, and the thrust coefficient 
decreases from this point onward. At the nozzle exit, the thrust coefficient 
is about 1.0, which indicates that the mixing gain has offset the decrease in 
thrust attributable to viscous losses. In turn, the thrust coefficient at 
the nozzle exit is 0.98 for the cold flow case, and thus the net gain in 
thrust due to mixing is about 2 percent. 

The final set of flow calculations considered here consists of Cases 5 
and 6, which correspond to lobe shape B-2 in Fig. 16. These cases represent 
cold flow conditions with equal fan and turbine static temperatures, but 
with different average inlet Mach numbers and Reynolds numbers. The ratio of 
fan to turbine inlet Mach number and velocity is 0.5 in each case. These 
flow conditions attain a large velocity difference between the fan and turbine 
streams as a result of mismatched total pressures and Mach numbers. Similar 
flow conditions are present in a related experimental study of flows in this 
same mixer geometry, being performed as part of this overall contract effort. 
These cases also serve to test the present calculation procedure under condi- 
tions of significant mismatch in fan and turbine total pressures and Mach 
numbers. The axial variation of Mach number and total pressure loss coefficient 
for Cases 5 and 6 are shown in Figs. 20 and 21. The total pressure losses in 
Fig. 21 are seen to be higher for the flow with higher initial Mach number, 
as is to be expected. These cases serve to demonstrate the capability of 
the present calculation procedure for treating flows with mismatched Mach 
number and total pressure. 
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SUMMARY AMD CONCLUSIONS 


A finite-difference method has been developed for detailed computation 
of three-dimensional subsonic turbulent flows in turbofan lobe mixers. 

The method is based on a decomposition of the velocity field into primary and 
secondary flow components which are determined by solution of equations 
governing primary SK>mentum, secondary vorticlty, thermal energy, and continuity. 
A two-equation turbulence model for turbulent kinetic energy and dissipation 
rate is used to construct the required Reynolds shear stress correlation 
coefficients. The governing equations are solved by a forward-marching solution 
procedure which corrects an ^ priori Inviscld potential flow solution for 
viscous and thermal effects, secondary flows, total pressure distortion and 
losses, internal flow blockage and pressure drop. 

Test calculations for a turbulent coaxial Jet flow were performed to 
verify that the turbulence model would perform satisfactorily for this 
relatively simple flow for which adequate experimental measurements are 
available. Several calculations for lobe mixer flows were performed for two 
geometries typical of current mixer design. These calculations included 
both hot and cold flow conditions, and both matched and mismatched Mach number 
and total pressure in the fan and turbine streams. The axial development of 
cross-sectional averages of Mach number, total pressure losses, and thrust 
were presented for these mixer flows, and these predictions are consistent 
with the behavior expected for this type of mixer flow. The results of this 
study serve to demonstrate the potential value of the present calciilatlon 
procedure to aid in the design and evaluation of turbofan lobe mixers. Further 
evaluation of the procedure will require careful comparison with sufficiently 
detailed and accurate experimental measurements. 
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TABLE I 

SUMMARY OF LOBE MIXER CALCULATIONS 


CASE P 

GEOMETRY 

T /T 
F' T 

F T 

Hj/Ht 

AVERAGE 
INLET M 

FAT^ STREAM 
Re 

LAMINAR OR 
TURBULENT 

1 

A 

0.5 

0.707 

1.0 

0.34 

1280 

L 

2 

A 

0.5 

0.707 

1.0 

0.34 

2.9x10^ 

T 

3 

B 

0.4 

.. . . ... 

' 

0.7 

1.107 

0.194 

2.6x10^ 

T 

U 

B 

1 

1.107 

1.107 

0.194 

2.6x10® 

T 

5 

B 

1 

0.5 

0.5 

0.0529 



4.75x10^ 

T 

6 

B 

1 

mm 

0.5 

0.171 

C 

1.54x10 

T 
























Figure 1. - Schematic of coordinate system for 
three-dimensional flow problem. 
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Figure 2 . 


. Typical lobe mixer geometry 


.„d potenti.1 flow coordinate «,ete«. 


Figure 


Casing 



3^ — Typical cross section at mixer exit surface. 
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Figure 4. - Family of lobe contours which may be specified. 
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Velocity deficit 



Figure 6. ■> Comparison of predicted and measured centerline velocity decay. 
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Computed 



Dimensionless axial distance, x/D 




100 


Figure 7. - Comparison of predicted and measured spread rate. 




Figure 8. - Comparison 


# Computed > x/D * 6.0 
■ Computed - x/D ■ 20 
▲ Computed - x/D ■ 37 

Data range of 

Forstall 6 Shapiro (Ref. 23) 


Note: Symbols Indicate 

grid points 


predicted and experimental velocity profiles. 





Figure 10. - Comparison of computed fine mesh velocity profile with 
assumed coarse mesh starting profile. 





Radius 


Casing 


of Calculation 



Figure 11- 


. computed coordinate system tor mlrer 


Geometry A. 






(not to scale) 


Figure 12. - Lobe shape at mixer exit surface 
for Geometry A. 
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Mass averaged total pressure loss coefficient. 



Dimensionless axial distance, x/L^ 


Figure 13. - Comparison of laminar and turbulent losses in mixer Geometry A. 



Thrust coefficient, T/T 



Figure 14. - Comparison of ideal thrust coefficient for laminar and turbulent 
flow In Geometry A. 




Figure 15. - Computed coordinate system for mixer Geometry B 
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Figure 17 



Area averaged Mach number 



Dimensionless axial distance, x/L^ 


Figure 18. - Comparison of axial Mach number distributions for hot and cold 
flow in mixer Geometry B. 
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Dimensionless axial distance, x/L 
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Figure 19. - Comparison of axial thrust development for hot and cold flow 
in mixer Geometry B. 
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Area averaged Mach number 


High initial M (Case 6) 
Low initial M (Case 5) 
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. - Comparison of average Mach number variations for flows 
with mismatched initial Mach numbers* 



Mass averaged total pressure 
loss coefficient, C„ 



Dimensionless axial distance, x/L 

» r 


Figure 21. - Comparison of losses for flows with mismatched initial Mach 
numbers and different average initial Mach number. 



APPEND IX-A 


THE TURBULENT PRODUCTION TERM 


The production term, P, in the turbulence kinetic energy and dissipation 
equations may be expressed as 


2(e„* + ej2'+e33‘') + e3*+e,/ + ej3*j 


where 




I ^ V ^ w_ 


dh. 


h| dx hjhj hjhj di 


(A. 2) 


•22 = 


; dv w ^ ^2 U ^^2 

h2 3 y hjhj d2 h|h2 


(A. 3) 


« _ _L + u V -^^3 

” hj az hjhj dx hgh, dy 


(A. 4) 


h 2 3y ^hs) H 7 ^2 1^2/ 


(A. 5) 


. ^ !!l Ll!L\ + liL i_ /J!L\ 

h3 dz Vh|/ h| dx \ h3/ 


(A. 6) 


_ i!L ^ (jl\ + hi JL / jl \ 

‘ hg dy \h, / h, dx ^h2/ 


(A. 7) 


The production term may be simplified considerably if it is first noted 
that in the present geometry the metric coefficients are not functions of z 
(the azimuthal coordinate). Furthermore, it may be assumed that the shear 
layers are thin, since the Reynolds number is high. Thus, if the flow is 
sheared in the y direction, v and 3/3y are assumed to be of order 6 and 
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1/6 » respectively. If the flow Is sheared in the z direction, w and 3/3z 
are assuned of order 6 and 1/6. The quantity 6 is taken to be representative 
of the shear layer thickness. Slinllarly, u, x and the netrlc coefficients 

are presuaed to be of the order of unity. If the viscosity, Is assumed 

2 ■*’ 
to be of the order 6 the relevant terms in F may be found by retaining only 

those terms which are of the order of unity In either of the two types of 

shear layer described above. As a result P is approximated as 



APPEND IX-B 


USER’S MANUAL - SRA MIXER CODE 


SUMMARY 

The use of a finite-difference procedure for computing turbulent flow in 
lobe mixers is described. The procedure is based on the solution of an approx- 
imate set of governing equations derived from primary and secondary flow 
concepts. Heat transfer is accounted for through introduction of an energy 
equation. Turbulence is modeled using a two-equation model. The computer code 
may be used either in a mode in which starting conditions are generated by 
automated procedures or in a mode where sufficiently smooth and detailed 
experimental data may be used to specify the initial conditions. The manual 
describes the analysis upon which the code is based, input required and output. 
The manual also contains brief descriptions of the subroutines which comprise 
the code as well as a list of principal Fortran symbols. A flow chart of the 
overall program algorithm is also Included. 
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INTRODUCTION 


The lobe mixer is a device currently being used for thrust augmentation 
on a variety of turbofan engines. The forced mixer is designed to mix primary 
(turbine) and secondary (fan) flow streams before they enter the exhaust nozzle, 
thus providing the nozzle with air having a more uniform energy which results 
in a more uniform velocity at the nozzle exit. It has long been known that 
by mixing the fan and turbine exhaust streams of a turbofan engine in this 
manner prior to expansion through the exhaust nozzle, a small but significant 
performance gain may be realized. The level of these gains depends on trade- 
offs between the degree of mixing of the two streams and the viscous losses 
incurred in the mixing process. 

To date, the performance of lobe mixers has been determined almost entirely 
through experiments. These experiments are sufficient to determine the relative 
merits of one mixer configuration as compared to another but they do not directly 
provide Information as to the reasons why one given mixer configuration performs 
better than another. The overall mixer flow field consists of two major flow 
regions; (1) the flow upstream of the lobe exit plane (within the lobes) and 
(2) the flow downstream of the lobe exit plane where the mixing actually takes 
place. The present study concentrates upon the actual mixing process downstream 
of the lobe exit plane, through development of a computational procedure 
applicable in this region and capable of predicting the three-dimensional 
mixing process in detail. 

Reference numbers and figure numbers mentioned in this appendix refer to 
the references and figures at the back of this appendix. 



ANALYSIS 


General 

The calculation procedure described here is an application of the 
approach developed by Briley and McDonald [1] to lobe silxer flows. The 
procedure is based on a deconposition of the velocity field into pr Inary and 
secondary flow velocities. Equations governing the streanwlse development of 
the primary and secondary flow velocity fields are solved by an efficient 
algorithm using both block and scalar ADI methods derived from the Douglas-* 
Gunn [2] splitting. Although the governing equations are solved by forward 
marching, elliptic effects due to curved geometries and area change are 
accounted for ^ priori through imposed pressure gradients determined from the 
potential flow in the geometry in question. Since the primary concern in the 
lobe mixer problem is thermal mixing of the fan and turbine stream to achieve 
thrust augmentation, an energy equation Is Introduced. Finally, turbulence 
is modeled using a two— equation k— e model [3]. Since a detailed explanation of 
the analysis and of Its application to the lobe mixer problem Is available 
elsewhere [1,4], only a brief discussion is given here. A flow chart of the 
calculation procedure Is shown In Fig. 1. 

The Governing Equations and Solution Procedure 

The solution procedure centers around the decomposition of the velocity 
Into primary and secondary flow components 

U-Up+U, (2.1) 

It Is assumed that in the orthogonal potential flow coordinate system used In 
the lobe mixer code, the primary flow velocity Is aligned with the potential 
flow streamlines (see Fig. 2). Similarly, the secondary velocity Is assumed 
to lie In the transverse coordinate surfaces which are normal to the stream- 
lines. The equation governing development of the primary velocity component 
is thus approximated as [4] 
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du du du 

dh. o dhq 


+ />(v,)uh3 


dy 


dm 

9 dh \ api dpjx) 

■*■ ^2^3 "air ^2^3 “d]; 


( 2 . 2 ) 


1 d 

■ (;i+/i,.|.)h,h 3 du 

h,hj a 

■ 1 

Re ay 

hg dy . 

Reh^ dz 

L dz J 


In Eq. (2.2) it is noted that the coordinate system is axisynnnetric and thus 

the metric coefficients are not functions of the azimuthal coordinate, z. 

Additionally, in Eq. (2.2) u represents the sole component of U , P_ is the 

P 1 

imposed pressure and P^(x) represents the mean viscous pressure drop which is 
determined from the integral mass flux condition 


j5fh,h,/)udydz 


» c 


(2.3) 


To account for thermal mixing between the fan and turbine streams, an energy 
equation is introduced in the form 

dE dE dE 



(2.4) 

with the assumption that the laminar and turbulent Prandtl numbers are equal 
to 1.0. 


The gas law is then introduced in the form 


d 

dm 



dPj dPy(x) 

dx dx 


(2.5) 


to relate the imposed and viscous pressure gradients to the other dependent 
variables. 
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At each axial station Eqs. (2.2) and (2.4) are solved initially using a 
block ADI scheme [5] with the convective coefficients lagged and with an 
assumed value for dP^(x)/dx. The density is then updated using Eq. (2.5) and 
the Integral mass flux is evaluated using Eq. (2.3). In general Eq. (2.3) 
will not be satisfied and Eq. (2.2) is thus solved assuming a nev value for 
dP^(x)/dx. This process is repeated using the standard secant method [6] to 
determine successive estimates for F^(x) until Eq. (2.3) is satisfied. At this 
point the priaury flov velocity component u, the total energy E and the density 
p have been determined at the new axial station and the secondary flow velocities 
must be updated. 

To determine the secondary flow velocity, is presumed to consist of 
solenoldal (rotational) and Irrotational components. The rotational components 
of are determined from solution of a streamwlse vortlcity equation [7] 
written in the following form [4] 


dC 
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Jii_/ 

d€ 

dpu 

dx 

^ dx 1 

0 i»2 V ** dy * dy / 
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h, 1 


dz 


2 dh, 

1 du h, I (■ dP 
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dp 

1 

(2.6) 



^ h, dz h,h, p 1 dy 

dz 

■ dV 

dz J 


P 

h| d 


h, d 

1 dlp.*p^)€ 1 

Re 

hjh, dy 

h,h,/) «y J 

h,* dz 

[p dz i 


2 2 2 

u ^8 vortlcity ( is assumed to lie in the transverse 


Here Q 

coordinate surfaces and is related to the velocity by 

dh.w. dh.v, 


I / fVs_\ 

« • h^h, \ dy ' dz I 


(2.7) 


Eq. (2.6) is solved using the Douglas-Gunn splitting of the Crank-Nicolson 
scheme. Once a solution for 6 Is obtained, a vector potential is determined 
from the equation 


' — f/ 

r \ 

1 1 ^ ^ » f/J 

1 \ a(h,^,)] 

h,h, ay H 

1 

' ay J dz [\ i 

£> / dz J " 


J --e (2.8) 


which is solved using iterative ADI. The rotational components of the secondary 

flow are then determined from Z as 

^s 
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(2.9a) 


_ I a(h,fe) 

V * />h,hj di 

* />h,hj dy 


(2.9b) 


At this point the resulting velocity field satisfies integral continuity but 
not local continuity. This is corrected through introduction of a scalar 
potential governed by the equation 

d(hjh3;)u) 


— f 

^h,h. 


d 

[ />h,^ 


dy 1 


dy J 

dz 

1 "j 

dz J 

i solution of 

Eq. (2.10) 

by 

scalar ADI 

J the ! 


dx 


( 2 . 10 ) 


nents v and w are found to be 

<J> <P 


IT w 


I ^<f>. 


and the secondary flow velocities are thus given by 

Vs • 


(2.11a) 

(2.11b) 


* w.,. + w 




(2.12a) 

(2.12b) 


The resulting velocity field 

u ■ ',Up + Uv, +T,W, (2.13) 

thus satisfies both differential and integral continuity. 
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The Turbulence Model 


The foregoing discussion was limited to mean flow equations and their 
solution. For turbulent flows, the turbulent viscosity which appears in 
Eqs. (2.2), (2.4) and (2.6) must be determined. The turbulent viscosity is 
determined using the two-equation k**c turbulence model as presented by Launder 
and Spalding [3]. In orthogonal coordinates and with the approximations 
required for forward marching solution [4], the equations governing the 
turbulence kinetic energy and dissipation are given by 

dk dk dk 

/^*'zh-Sr * A>v,h,hs-5^ + /»w,h,hg-5j- 


' -i-f 

h,hs 


1 — 1 

Re dy 1 



[ dyJ 


® r/ Mt \ 


C* € 

- hjhgh^Cg^-jj— + Cj-jj-p 

where k is the turbulence kinetic energy, e is the turbulence dissipation and 
P represents the turbulence production given as 

r/* . /I dw? . /I du\^ , / I dv\^i 


f/i ouv*" /I dwf /I du\ . /I dvri 

^Tl(h 2 dy/ 'h 2 dy' (h^ dz) (hsdz/J 


(2.16) 


Eqs. (2.14) and (2.15) are solved independent of the fluids equations using 
a linearized block implicit (LBl) scheme. The turbulent viscosity is then 
determined from k and c by the relationship 

R. (2.17) 
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The constants appearing in £qs. (2.14) » (2.15) and (2.17) are assigned their 

usual values [3] of; c_ = 1.44, c« = 1.92, c = 0.09, cj = 1.0 and a = 1.3. 

X z y £ 

It should be pointed out that the production term in the k and c equations 
may be expressed computationally in several ways which may be selected by the 
user. The first option available to the user is to have the production terms 
expressed as in Eq. (2.16) with the turbulent viscosity lagged. In a second 
option, the production term in the dissipation equation is altered by substi- 
tuting the relationship of Eq. (2.17) for y^. This simplifies the term 
representing production of dissipation and the resulting expression appears to 
be computationally advantageous. In the third option, the substitution for 
y^ in the production term is made in both the k and c equations. Computationally, 
the second option appears most attractive and is presently recommended. 

If all the available two-equation turbulence model options are exhausted 
and a successful calculation is still not obtained due to turbulence model 
failures, the user may select a simpler wake turbulence model. The use of 
the wake model should provide a means of completing a given calculation under 
these circumstances 


The Coordinate System and Specification of the Imposed Pressure Field 


The coordinate system used in the present MIXER code is generated using 
the ADD code [8]. The ADD code generates a two-dimensional potential flow 
coordinate system, an example of which is shown in Fig. 2. The axisymmetrlc 
coordinate system for the lobe mixer calculation is then obtained by rotating 
this two-dimensional coordinate system about a centerline. Although the 
resulting coordinate system is suitable for the viscous flow calculation [4], 
the two-dimensional pressure field is not a suitable approximation to the 
axlsymmetric pressure field. An a priori axisymmetrlc potential flow calcula- 
tion is performed to obtain a suitable pressure field. The axlsymmetric 
potential is obtained from the solution of the equation 


h,t,3 a f 

hgh, 1 d<f> h,hj 
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h, dx^ [ 

h, J dx hj 

dy^ dy 
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\ ay 


O (2.18) 


using an iterative ADI scheme. The resulting invlscid axisymmetrlc velocity 
field is then used to determine Incompressible axisymmetrlc pressure 
gradients 
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(2.19) 


( UfUi \ - I 


The streanvise pressure gradient is then scaled for compressibility using the 
local Mach number, M. , which is determined as part of the viscous flow 


calculation, to give the 


imposed pressure gradients 

ax * dx 


( 2 . 20 ) 


The maximum value of used in £q. (2.20) is limited to 0.7 in the computer 
code to avoid a singularity at * 1.0, although local Mach numbers greater 
than 1.0 may appear in the flow field [4]. 
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SPECIFICATION OF INITIAL CONDITIONS 


The initial conditions for the lobe MIXER code are specified as input, 
with three possible starting procedures. The IC0AX=1 option is used for 
simple two-dimensional test cases. Under the IC0AX=1 option a confined 
coaxial jet Initial profile is constructed assuming a cosine velocity 
profile for the shear layer between the primary and secondary streams (not 
to be confused with the primary and secondary flow) as shown in Fig. 3. The 
velocity and static temperature ratio between the two streams is specified 
as input by the user. A lobe of constant radius is specified, and the cosine 
shear layer velocity profile is centered about the lobe radius. The thickness 
of the shear layer is input as the variable DELMIX. Boundary layers may be 
specified at and/or (see Fig. 3) using the IHUB and/or ISHR options. 

Turbulence quantities are initialized by specification of a length scale which 
is assumed to be constant across the shear layer but which varies with distance 
from a wall if boundary layers are specified [4]. Additionally, free stream 
values of the turbulence kinetic energy are specified by the user. The user 
also has the option of initializing turbulence quantities by specifying constant 
turbulence kinetic energy throughout the flow field by using the IDIS option. 

A second method of starting the calculation procedure is available to allow 
the user to perform calculations using automated starting profiles generated by 
the code itself. This starting procedure is used unless either the ICOAX option 
or IDATA option (to be discussed subsequently) is requested. In this mode, 
the user must specify a lobe shape, in terms of discrete data points, as 
indicated in Fig. 4, starting with the minimum r coordinate value. The lobe 
coordinates must be specified such that the lobe enters the computational 
domain through boundaries 1 or 4 and exits the domain through boundaries 2 or 
3 as shown in Fig. 4. The user has the option of specifying boundary layers on 
the lobe, hub and shroud surfaces, or any combination thereof through use of 
the ILOBE, IHUB and ISHR options. The primary (turbine) and secondary (fan) 
streams are indicated in Fig. 4. The user may specify nominal velocity and 
static temperature ratios between the streams. Initial velocity profiles are 
constructed by assigning the appropriate values for velocity to each scream. 

The resulting uniform profiles are then corrected for radial pressure gradients 
present at the initial plane as determined from the potential flow calculation 
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and are also scaled by boundary layer profiles. Turbulence quantities are 
again initialized through specification of a length scale and free stream 
tturbulence intensities. However, the user may specify whether the length 
scale in the region near the lobe varies with distance from the lobe or Is 
constant across the lobe wake. Near the hub and shroud, the length scale will 
vary %rlth distance from the wall If botmdary layers are specified. 

The third starting procedure is activated through the IDATA option. If 
IDATA is greater than zero, the program accepts card input to specify the 
velocity and temperature field at the initial plane. The data must be supplied 
for each grid point in the computational domain. A lobe shape is specified as 
before, and DKT.MIX is input as a number representative of the boundary layer 
thickness on the lobe. Boundary layers on the hub and/or shroud may be 
considered using the IHUB and ISHR options. If these options are elected, 
DELHUB and DELSHR must be input as representative values of the boundary layer 
thickness on the hub and shroud, respectively. Turbulence quantities are 
initialized as in the previously described starting method. 
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DESCRIPTION OF INPUT 


Input to the MIXER code is accomplished by a combination of card input 
and information which is stored on auxiliary files. Metric information 
generated by the ADD code [8] and restart data sets are stored on separate 
auxiliary files. Program control is performed through card input, which is 
divided into three categories: 

1) Plot file information for NASA-Lewis plotting package, 

2) NAMELIST input, and 

3) Experimental flow field input. 

Format of the required input is as follows; 

Plot File Input 


CARD 

if 

COLUMNS 

VARIABLE 

DESCRIPTION 

Card 

1 

Col. 1-32 

TITLE (I) 

Plot title, format 5A6,A2. (May be 
left blank) . 

Card 

2 

Col. 1-2 

ISYN 

Twice the number of lobes in a 360° 
cross section of the mixer geometry. 



Col. 3-12 

SYSTEM 

Indicates the type of coordinate system 


for plot routines. Format FIO.O, 
Use SYSTEM=2. 


Namelist Input 


Namelist 

$IN1 

Description 

RG 

Gas constant, ft lb^/(slug °R) , default 1716.3 (air). 

CP 

Constant pressure specific heat, ft lb -/(slug °R) , 
default is 6012.384 (air). 

TZERO 

Reference temperature, assumed to be nominal primary 
stream temperature in degrees Rankine. 

UZERO 

Reference velocity, assumed to be nominal primary 
stream velocity, ft/sec. 

PZEROl 

Reference static pressure, assumed to be nominal 
pressure at initial plane, Ib^/ft^, default is 2116.8. 

YZERO 

Reference length, ft. 
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Namelist 


Description 


VISCl 

$END 

$ZN2 

lAXI 


ICOMP 


IDATA 


XEI'ITR 


YS(I,J) 


Reference viscosity, slugs/ (ft sec) . 


■1, cylindrical coordinates. 

•2, rotated orthogonal coordinates (from ADD code) 
(default). 

>■-1, stop after input dump. 

«0, normal run (default) . 

*1, stop after printing initial plane geometry. 

«2, stop after printing initial flov field. 

*0, automated start (default) . 

■1, read initial velocity and temperature field data, 
use wall length scale near lobe. 

■2, read Initial velocity and temperature field data, 
use wake length scale near lobe. 

Dimensionless distance from initial ADD code potential 
line to initial plane of lobe mixer calculation. 

Used only for lAXl-2. Determined from ADD code as 
XENTR" ( JSTE Pyg^r pi~l) *DS where JSTEPxentR number 

of the ADD code potential surface representing the 
MIXER code initial plane. See Fig. 2. 

Dimensionless transverse limits of computations domain. 

YS(l,l)-rjjjjj 

YS(2,l)»rj^ 

YS(l,2)»0jjjjj (radians) 

YS(2,2)»0jj^ (radians) 

See Fig. 2 for definition of YS(I,1) for mixer geometry. 
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Namelist 

NE(I) 


IGRID(I) 


EPS (I) 


XCTR(I) 

NS 

AP 


ILAP 


LX 


Description 

N£(l) * number of radial grid points. 

For computations using ADD code geometry N£(l) must 
be equal to the number of ADD code streamlines. 

N£(2) B number of azimuthal grid points. 

Maximum values for N£(l) are 20. 

IGRID=0, no grid stretching (default). 


= 1 , 

- 2 . 
or 0 


Roberts stretching at or 

Roberts stretching at both and 

MIN ®MAX' 

hyperbolic sine stretch about XCTR(I). 


1=1 for radial direction, 1=2 for azimuthal direction. 

Used for IGRID(I)?^0 to control the grid stretch. 

If unequally spaced streamlines are used in the ADD 
code geometry, the value of EPS(l) must be specified 
from the ADD code. 

For Roberts stretching 0.0<EPS(I)<1.0. 

For hyperbolic sine 1.0<EPS(I)<5.0. 

Used only with IGRID(I)=4 to specify the r(I=l) or 
6(1=2) location about which points are concentrated. 


Number of axial stations to be computed. 

NS=1 is initial plane. Maximum value is 50. 


If AP is input .GT.O a geometric streamwise grid is 
set up with Ax(J)/Ax(J-l) = AP. 

If AP is input .LT.O the user must specify the axial 
mesh. AP has a default of 1.05 but the AP.LT.O option 
is recommended. 


ILAP=1 Indicates that an axlsymmetric potential solu- 
tion is to be computed on the ADD code grid and the 
resulting pressure field used in the viscous flow 
calculation. ILAP has a default of 1. If lAXI is 
input as 1, ILAP is set to 0. 

LX is the number of potential lines generated by the 
ADD code. LX is used with ILAF=1. The maximum 
allowable value of LX is 50. 
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Hamellst 

LY 


IPCOR 


LSEC 


ILAM 


Description 

LY is the number of streamlines generated by the 
ADD code. LY is used with ILAF*1. The 
allowable value of LY Is 20. 

IPCOR is used with 1LAF>1 to identify the potential 
surface whose average dimensionless velocity is 1.0 
in the axisymmetric potential calculation. IPCOR 
should be input as IPCOR* JST£P__j_ in general. 

See Fig. 2. 

LSEC is a secondary flow option. 

LSEC*1, solve streamwise vortlclty equation. 

LSEC*0, do not solve streamwise vortlclty. 

LSEC has a default value of 1. 

ILAM*0 for turbulent flow, default. 

ILAM"1 for laminar flow. 

ILAM»-1 wake turbulence model, requires 
specification of freestream turbulence for 
stream 1 and a length scale. 


The default value is 0. Only laminar flow may be 
considered using the automated starting routine 
without boundary layers. 

ICOAX lCOAX-0, default. 

IC0AX*1 activates the coaxial Jet starting option. 

X(JX) X is an array which contains NS values of the stream- 

wise mesh coordinate. 

If AP is input .GT.O only X(l) and X(2) need be input. 

If AP is input .LT.O NS values of X must be input. 

If IAXI*1, X is a physical but dimensionless distance. 

If IAXI*2 and ADD code geometry is used, X represents 
the dimensionless computational distance from the 
Initial plane. 

The maximum allowable value of X is given as 
(JSTEPmaX"*^STEPxenTR>*®S/USCALE where USCALE-YS(2,1)/ 
(YS(2,1)-YS(1,1)) . The maximum value of X should be 
slightly less than this value to avoid marching out of 
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Namelist 


Description 


the ADD code geometry which will result in an error 
termination. See Fig. 2. 


IKECUP 


IDIF 


IDIFBC 


IDIS 


$END 

$IN3 

IRSTIN 

IRSTOT 


JRSTIN 

JRSTOT 


*0, lagged In turbulence model equations production 
terms. 

*1, treated Implicitly in turbulence dissipation 
equation (default) . 

■2, treated Implicitly In both k & c equations. 

*=0, central difference radial convective terms in 
turbulence model equations. 

»1, use two-point one-sided differencing of radial 
convective terms in turbulence model equations (default) . 

=0, use three-point one-sided normal gradient boundary 
condition weights. 

=1, use two-point one-sided normal gradient boundary 
condition weights (default) . 

=1, specify uniform turbulence kinetic energy, at the 
turbine stream value throughout the initial flow field, 

=0, initial turbulence kinetic energy assumes 
turbulence production equal to dissipation. 

IDIS=0 by default. IDIS*1 should be used only with 
the constant length scale options th^s if IDIS=1 and 
ILOBE is input as 1, ILOBE will be set to 2. If IDATA 
is input as 1, it will be set to 2. 


The number of the axial station to be read in for a 
restart. No restart if IRSTIN=0 (default). 

Increment for saving restart information. Restart 
files will be written every IRSTOT steps from starting 
value of JX, the axial station counter. Default 
assumes no restarts written. 

File number from which restart information is to be 
read. Default file number is 11. 

File number on to which restart information 
is to be written. Default file number is 11. 
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Hamellst 

NFILE 

HSAVED 

IPLOT 

$END 

$IN4 

USTRM2 

TSTRM2 

TUSTRl 

TUSTR2 

ALEN 

AMEXIT 

DELMIX 


Description 

Sequence nua^er In file JRSTIN of the desired restart 
infomatlon. 

The ntsnber of restart data blocks saved on file 
JR5T0T. 

If IFLOT-l a plot file must be assigned in the run 
stream and Information will be written In this file 
for subsequent use in plotting. IPLOT>0 give no 
plot file. Default is IFLOT-0. 


Ratio of secondary stream (fan) velocity to primary 
stream (turbine) velocity. For matched a verage i nlet 
Mach numbers and total pressures USTRM2* ^TSTRM2 . 
Default value Is USTRM2*.707. 

Ratio of secondary stream temperature to primary stream 
temperature. Default Is TSTRM2*.5. 

Primary stream turbulence Intensity. 

Tuj^- u*^/uj where u^»UZERO. 

Secondary stream turbulence Intensity. 

T„ 2 “ u*2/u 2 where u^-UZEROxUSTRMZ . 

If TDSTR2 is Input as 0.0 the free stream turbulence 
In both streams (k-3T^2/2) is assumed equal to 3Tu^/2. 

Dimensionless free stream length scale l^^^/YZERO. 

To be estimated by the user. As a guide. If IDATA or 
ICOAX is Input .GX.O ALEN should be on the order of 
O.lxDELMlX. If the automated MIXER starting routine 
Is used ALEN should be on the order of the displacement 
thickness specified. 

AMEXIT is an assumed exit plane Mach number which Is 
used to specify the Isentroplc exit plane pressure for 
the thrust calculation. Default value Is AMEXIT»0.9. 

DELMIX is used in both the IDATA and ICOAX options. 

If IDATA is input .GT.O DELMIX must be input as a 
dimensionless boundary layer thickness, d/YZERO, 
presumed to be representative of the boundary layer 
on the lobe surfaces. If ICOAX is input .GT.O 
DELMIX is the dimensionless thickness of the cosine 
shear layer. See Fig. 3. 
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Namelist 

DELmiB 

DELSHR 

$END 

$IN5 

ILOBE 


IHUB 


ISHR 

BLHUB(I,J) 


Description 

DELHUB is input with IDATA .GT.O as a dimensionless 
boundary layer thickness, 6/YZERO, representative of 
the hub boundary layer. 

Similar to DELHUB but pertaining to the shroud 
boundary layer. 


ILOBE is a flag used in the automated starting 
routine. 

1LOBE=0 indicates that no boundary layers are to be 
set up on the lobe surfaces. 

IL0BE=1 indicates that lobe boundary layers are to be 
set up and a wall-type length scale based on distance 
from the lobe is used to initialize turbulence 
quantities. 

ILOBE-2 is similar to ^1 but the length scale is held 
constant across the lobe wake. 

If ILOBE is input .GT.O the IL0BE=2 option is 
recommended. If IDIS is input as =1 and automated 
lobe boundary layers are to be constructed, ILOBE will 
be reset to 2. 

IHUB=0 indicates that hub boundary layers will not be 
present in the initial profile, regardless of other 
starting options. If IDATA is input .GT.O and IHUB 
is input .GT.O DELHUB must be input. If the IDATA 
option is not used, IHUB must be input equal to the 
number of points input to describe the hub boundary 
layer in the array BLHUB. 

ISHR is similar to IHUB except it pertains to the 
shroud boundary layer. 

If hub and/or shroud boundary layers are to be 
constructed using the automated starting routine 
(IDATA^O) the boundary layer displacement thickness 
and shape factor must be input as functions of 6 along 
the hub and/or shroud surfaces. This is accomplished 
using the array BLHUB. 
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Hamellst 


Daacrlptlon 


BLHUB(Z,1) contains the values of 0 (in radians) 
where hub boundary layer paraaieters are specified. 

The 6 values are arbitrary but aust span the computa- 
tional domain. The number of 6 values specified snist 
equal the input value of IHUB. A maximum of 20 values 
may be specified. 

BLHUB(I,2)*6*/YZERO» the dimensionless displacement 
thickness of the hub boundary layer at each specified 
6 location. 

BLHUB(I»2)*H, the shape factor at each 6 location. 

BLHUB(I,4) la the shroud equivalent of BLHUB(I,1). 

ISHR values of 6 must be loaded in BLHUB(I,4). 

BLHUB(I,5) and BLH0B(I,6) are similar to BLHUB(1,2) 
and BLHUB(I,3) except they pertain to the shroud 
boundary layer parameters. See Fig. 4. 

NLOBE NLOBE is the number of data points input to describe 

the lobe shape and lobe boundary layers. 

YL0BE(I,J) YLOBE is an array which contains the data needed 

to specify the lobe shape and lobe boundary layers. 
NLOBE values of each parameter in YLOBE must be read 
in. 


YLOBE (1,1) contains NLOBE dimensionless radial 
coordinates, R/YZERO, of the points used to describe 
the lobe shape, starting at the minimum radius. 

YLOBE (I, 2) contains the corresponding 6 coordinates 
in radians. 

YLOBE (I, 3) contains NLOBE values of the dimensionless 
displacement thickness, 6*/YZ£R0, at each coordinate 
point, for the primary stream lobe boundary layer. 

YLOBE (I, 4) contains NLOBE values of the displacement 
thickness for the secondary stream. 

YL0BE(I,5) contains the distribution of the shape 
factors, H, at each coordinate point for the primary 
stream. 

Finally, YL0BE(I,6) contains the shape factor distri- 
bution for the secondary stream. 
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Namelist 


IWAKE 


$END 

$IN6 

NBCON(I,J,K) 


Description 

YL0BE(I,3) through YL0BE(I,6) need not be specified 
if the IDATA or ICOAX option is specified, or if 
ILOBE*=0. YLOBE(I,p and YL0BE(I,2), the r-0 coordinates 
of the lobe, must be specified for all starting options. 
See Fig. 4. 

If IWAK£>=1 the slots for shape factor in BLHUB and 
YLOBE need not be input as an assumed value of the wake 
function in the boundary layer profile is assumed. 


NBCON is a three-dimensional array which allows the 
user to select the desired boundary conditions for the 
streamwise velocity, the total enthalpy, the turbulence 
kinetic energy and the turbulence dissipation. It 
should be recalled that the momentum and energy equations 
are solved as a system of equations and the turbulence 
equations as a second system of equations. The third 
subscript, K, indicates which set of equations boundary 
conditions are being input for. 


K=1 for the momentum and energy equations. K=2 for the 
turbulence model equations. With R^^l, I«1 Indicates 
a boundary condition for u is to be specified. K=1 
and 1=2 is total enthalpy. K*=2, 1=1 is turbulence 
energy and K=2,I=2 is turbulence dissipation. 

The value of J ranges from 1 to 4 indicating the 
computational boundary considered. J=1 for the 

boundary at RmIN» ^MAX* ®MIN 

J=4 for 


NBCON(I, J,K)=0 specifies zero function value boundary 
condition. 


NBC0N(I,J,K)-1 

NBC0N(I,J,K)=2 

NBC0N(I,J,K)-4 

conditions. 


specifies zero normal gradient, 
specifies second derivative zero, 
specifies wall function boundary 


Default values are as follows: 

For streamwise velocity and energy symmetry (zero 
gradient) at and and zero second derivative 

at Rmin ^nd turbulence variables, normal 

derivatives are specified zero on all surfaces. 
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Nanellst 


Description 


XCBC 


MBCONC(l.J.R) 


$END 

$IN7 

IPRN(I) 


If a calculation Is being performed with hub/or shroud 
boundary layers (using wall functions) the velocity 
boundary conditions at 

specified as 4. Heat transfer at”^e wall Is not 
considered In the MIXER code at present thus NBC0N*1 
should be specified at %IAX total 

enthalpy. Default values should be used for the 
turbulence energy (which is consistent with wall 
functions) and MBC0N«4 should be specified at Bmxn 
R j,j^ for the turbulent dissipation. If hub and/or 
shroud boundary layers are not specified, the default 
values at and/or will suffice. 


XCBC is used In MIXER calculations to change boundary 
conditions as the flow moves down stream off the hub, 
irtiere wall function boundary are applied, to a center 
line, where symmetry conditions are applied. 

The user must specify the value of XCBC by examination 
of the geometry from the ADD code. XCBC Is computed 
as XCBC* ( JSTEPj^Qg^ ” JSTEPj^jj^j^) ^DS where 
the number of the ADD code potential line where the 
plug degenerates to a center line and JSTEP yirTi rnt 
number of the potential line representing the Initial 
MIXER calculation plane. See Fig. 1. The default 
value causes no changes. 


NBCONC has the same function as NBCON but is used to 
specify which boundary conditions are to be changed 
for X>XCBC and what they are to be. Default value 
will cause no change In boundary conditions regardless 
of the value of XCBC. For a calculation with hub 
boundary layer, NBCONC Is used to change from wall 
functions to synmetry conditions as the hub degenerated 
to a center line. 


IPRN Is an array which selects variables to be printed 
at each axial station. If IFRN(I)*1 variable “I" will 
be printed. If IPRN(I)*0 variable "I" will not be 
printed. The variables are numbered as follows: 
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Namelist 


Description 


Ib Variable 

1 u 

2 V 

3 w 

4 P 

5 E 

6 Pj 

7 K 

8 e 

9 e 

10 ij/ 

11 M 

12 static temperature ratio T/T^ 

13 total temperature T°/T° 

14 pressure coefficient 

15 static pressure ratio P/P^. 

16 total pressure ratio P*^/P^ 

17 swirl angle, degrees 

18 u'v’ 

19 

20 v*w* 

ICON ICON controls printer contour plots of the printed 

flow field variable. 

IC0N=»1 produces printer contour plots. 

ICON=0, no contour plots. 

$£ND 


84 



Experlaental Flow Field Input 


The follcwlng cards are needed only 1£ starting profiles are to be read 
in (U)ATA«1 or 2). Each card contains velocity components u, v, w in feet 
per second and static temperature, T, in degrees Rankine at each computational 
grid point. One card per grid point is required, thus a total of NE(1)*NE(2) 
data cards are needed. Format 4F10.0. The data is read for the grid point 
numbering shown in Fig. 5. 

CARDS # $ENDfl to $END^NE(1)*NE(2) 

Col. 1>10 u - velocity component in the computational streamwlse 
direction, ft/sec. 

11-20 V - velocity component in the computational radial direction, 
ft/sec. 

21-30 V - velocity component in the computational azimuthal direction, 
ft/sec. 

31-40 T - static temperature, ^R. 


Sample Run Stream 

Listed In Fig. 6 is a sample UNIVAC 1110 run stream for a lobe mixer 
calculation. The run stream includes all file assignments required excluding a 
plot file, using default values of JRSTOT and JRSTIN. The input data will 
reproduce the Case 6 results presented in [4] and will duplicate the sample 
output contained in this manual. 

The first card in the run stream assigns file 9. File 9 is a temporary 
file which is used to store out of core flow field variables. File FDATA5(2) 
is the file which contains geometric information. It is created using the ADD 
code [8] and must be a FASTRAN format file. The USE card which follows the 
assignment specifies the Internal unit number used in the MIXER code for the 
geometry file. PEPR3 is the restart file and snist be cataloged as a word 
addressable file before the run. The USE card which follows specifies the \mlt 
number which is used in the MIXER code for the restart file. Note that this 
must be the saaie as the value of JRSTOT. If JRSTOTi^ JRSTIN an additional file 
aust be assigned with a USE card specifying the unit number as JRSTIN. Files 12 
through 14 are temporary files. The remainder of the input has been described 
previously. 
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DESCRIPTION OF THE OUTPUT 


The output of the MIXER code can be divided into two categories. The 
first Includes output associated with initiation or restart of a calculation. 
The second category of output is controlled by the user and consists of 
stnmnary tables of Integral properties, skin friction coefficients, and flow 
field printouts at each marching station. Sample output, typical of that 
obtained when initiating a calculation using the automated starting procedure 
is shown in Figs. 7-13. 

The output for any case begins with the printed message ***INPUT NAMELIST 
DUMP*** which is then followed by a dump of namelists INI through IN7. This 
provides the user with a record of the namelist input used to run a given 
case. It should be noted that the values of XENTR and XCBC printed here may 
differ from the input values since they have been divided by the metric scale 
factor for consistency with the MIXER coordinates. After the namelist dump, 
the message ****SUMMARY OF INPUT DATA**** is printed. This summary data is 
shown in Fig. 8 and includes the R-0 coordinates which describe the lobe 
shape, input which describes the mixing duct including the computed value of 
the metric scale factor, the primary or turbine stream, reference conditions, 
the velocity, temperature and total pressure ratios between streams, free 
stream turbulence intensities and length scale, and a list of computational 
options either specified or set by default. Following the list of computational 
options. Information pertaining to the computational coordinates is printed. 

This information Includes NS values of the axial mesh distribution, the 
transverse grid limits and grid stretching Information. The message 
**TRANSVERSE GRID AND DIFFERENCE WEIGHTS** is then printed followed by a table 
as shown in Fig. 9. In this table, the first column contains the grid point 
Indices for the y or radial direction, followed by the z or azimuthal direction. 
The second column lists the computational location of each grid point. The 
third through fifth columns list the first-derivative difference weights and 
columns 6 through 8 list second-derivative weights. If the case is a restart, 
a restart message will be printed at this point and the calculation will 
proceed. If the case is an initial calculation and if the IAXI*=2 option with 
the ILAP=1 option is specified, the message **POTENTIAL FLOW SOLUTION** will 
appear next, followed by an ADI convergence message. The next LX columns of 
information printed is the axlsymmetric potential flow solution followed by 
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XJC columns of informatloa which represent the incompressible potential flow 
pressure coefficients. The potential flow solution will not be printed if 
IAXI-1. If IAXI»2 and ILAP«0, or if 1R5TIN>0. The potential flow printout is 
shown in Fig. 10. 

If IAXI>2, metric information will follow the potential flow solution (if 
it has been printed) or restart message. This metric information is shown in 
Fig. 11. JSTEF indicates the number of the potential surface from the ADD 
code Just upstream of the current axial location. DSTEF is the MIXER equivalent 
of the ADD code DS divided by the metric scale factor. SQl is the value of the 
axial coordinate at the JSTEF+1 potential surface and SQ2 is the axial coordinate 
value at the JSTEF potential surface. SQ12 is the value of the axial coordinate 
of the transverse coiiq>utatlonal plane for which the metric information is 
needed. The variable FRACT is the fraction of the distance between SQl and SQ2 
of SQ12, FRACT-<SQ12-SQ2)/(SQ2-SQ1). In the table which follows, TPHYS and 
XPHYS are dimensionless radial and axial coordinates in physical space, 
respectively. HI, H2 and H3 are the metric coefficients, DHl/DY, etc., are 
derivatives of the metrics and PRESS is 1/2 of the inviscld pressure coefficient. 
If this is the first geometry printout of a new case or restart, two tables 
of this metric information are printed. The first corresponds to the initial 
or restart location and the second corresponds to the following axial station. 

If the run is a new case having IDATA>0, for which Initial profile data is 
read in, the message ****UNFR0CESSED INPUT DATA**** will be printed followed 
by flow field prints of normalized specified input profiles. If automated 
starting procedures are used, this printout is omitted. The next item printed 
is the message ****AR£A INTEGRALS**** followed by a summary table of area>- 
weighted or averaged quantities. This summary table is shown in Fig. 12 and may 
be repeated several times to trace the history of the iteration of the mean 
pressure drop. The last table printed contains the values associated with the 
converged solution. The table contains, in dimensionless form, the following: 

AREA - cross->sectlonal flow area. 

MASS FLUX - self explanatory. 

VELOCITY FLUX - the area-averaged streamwlse velocity. 

AVE BD - the area-averaged value of pu. 

AREA RATIO - the ratio of the local area to the inlet area. 
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MASS ERROR - the error in the mass flux from the mean pressure drop 
iteration. 

MASS AVE CPT - a mass-averaged total pressure loss (cf. [4]). 

IDEAL THRUST COEF — a ratio of the local value of the thrust obtainable 
to the value at the mixer inlet (cf. [4j). 

2 2 

THRUS - the dimensionless thrust, T/(p u y ). 

r r o 

THARAT - an area ratio of the required exit area to the exit area required 
at the initial plane. This also implies losses as it rises 
above 1.0. 

THAREA - the exit area used in the thrust calculation (cf. [4]). 

AVE MACH NO - the area-averaged value of the Mach number. 

AVE PT RATIO - the area-averaged total pressure ratio. 

AVE T TOTAL RATIO - the area-averaged total temperature ratio. 

AVE CP — the area-averaged pressure coefficient. 

BLOCKAGE - a local reference quantity. 

Following this summary table, the message ****PR0CESSED STARTING PROFILES**** 
is printed as shown in Fig. 13 followed by flow field prints of the initial 
conditions after some preprocessing to be compatible with continuity. From 
this point onward, the form of the program output is independent of the 
mode of operation. With the exception of the first computational station, 
a table of geometric information is printed if IAXI=2. This is followed by 
the area integral tables. Following the area integrals, an ADI convergence 
message is printed indicating the convergence of the vector potential equation, 

Eq. (2.8). A second ADI convergence message is printed indicating convergence 
of the scalar potential equation, Eq. (2.10). A tabulation of skin friction 
coefficients is then printed if hub and/or shroud boundary layers are specified. 
Following the skin friction coefficients, the detailed flow field variables 
specified by the user are printed. The columns labeled R and Z (see Fig. 13) 
contain the dimensionless radial and axial location of each grid point in 
physical space. There is no azimuthal (6) variation of these coordinates 
due to the axisyrametric coordinate system. After the flow variables are 
printed, the output sequence is repeated for subsequent axial steps. 
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TEMPORARY Ain> PERMANENT STORAGE REQUIREMENTS 

With the MIXER code dimensioned for e 20x20 grid in the transverse 
plane and 50 axial stations, the core requirement on a DNIVAC 1110 computer 
is about 68000 words. Four temporary file assignments are required. These 
files are assigned as word-addressable files on the DNIVAC 1110 as follows: 

Unit 9 50,000 words 

Unit 12 5,000 words 

Unit 13 5,000 words 

Unit 14 5,000 words 

If a plot file is to be written, a permanent file must be created prior to the 
run. This file should be FASTRAN-formatted and must be given a local file name 
of Unit 8. The actual size of the file will depend on the computational grid 
used but approximately 14xNE(l)*NE(2)+5 words are written on this file for 
each axial step. If restarts are generated, a permanent word-addressable file 
must be created and a local file name consistent with JRSTOT must be used (unit 
11 by default). Again, the actual length of this file will depend on the 
number of restarts written but approximately 22000 words per restart will be 
written with the deck dimensioned as stated above. 


ERROR CONDITIONS 


Certain types of error conditions may cause program termination. Those 
conditions which are known to represent possible modes of error termination are 
discussed below. 

If both ICOAX and IDATA are greater than zero* the computer code will 
produce the message ***INPUT FOR IDATA AND ICOAX IS INCONSISTENT*** and the 
run will be terminated. The user should correct the input. 

Several error terminations may occur in subroutine GEOTRB. The first of 
these errors terminates the run with a RETURN 0 CALL FROM GEOTRB and is 
associated with the Inability of the computer program to find the position in 
the geometry file which corresponds to the current axial location to be 
computed. Error termination also occurs if the value of the axial coordinate 
is greater than the extent of the geometry generated by the ADD code. Both 
of these errors occur only for IAXI=2 and are caused by improper specification 
of the axial grid. The user should examine the grid being used and correct 
it as necessary. Error termination may also occur in GEOTRB if the geometry 
file is improperly assigned. 

Although no error termination occurs, it is possible that the ADI routine 
used to solve the vector potential equation or the scalar potential equation, 
Eqs. (2.8) and (2.10), may fail to satisfy certain convergence criteria. 

When this happens, it usually indicates a failure to meet the convergence 
requirement within the allotted number of iteration, which is set at 50, rather 
than indicating divergence of the solution. Should this situation occur, the 
message ***ADI FAILS TO CONVERGE*** is printed followed by the message 
ITERS, PHIMAX,RHSMAX,EPS1-2=XX .XXX .XXX .XXX .XXX 

where the X's indicate numbers. The first number is the maximum number of 
Iterations. The following two numbers represent maximum normalized errors 
in the solution. If the first of these two numbers is less than 1.0 and the 
second number is not greater than about 10., the message may be ignored although 
if repeated failure occurs, the secondary flow velocities should be examined 
carefully. If both of these numbers are significantly greater than 1.0, it 
may indicate divergence of the solution, although this has not been encountered 
as of this writing. Should such a failure occur, the user should carefully 
examine all input to determine if a problem exists elsewhere. The ADI messages 
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appear Just after the area Integral suimiiary tables. The first ADI message 
corresponds to the vector potential; the second corresponds to the scalar 
potential. 

Another error message is generated by the code (but does not terminate the 
run) if the iteration for the mean viscous pressure drop falls to converge within 
5 iterations. If this occurs, the message WARNIHG^MASS FLUX ITERATION DIO NOT 
CONVERGE~CHECR RESULTS CAREFULLY appears. If this message appears, the user 
should check the mass error printed in the area integral summary tables and 
decide whether the solution is acceptable. 

If the flow chokes upstream of the last exit plane, the warning message 
discussed above may appear at the axial station before the flow chokes. At 
the axial station where the flow actually chokes the calculation procedure may 
break down completely and an error termination may occur in the square-root 
routine during the iteration for the mean pressure drop. Under these conditions, 
the user should adjust the inlet mass flux to allow the calculation to proceed 
to the exit plane without choking. 

Program execution will also be terminated on a restart case if the input 
value of IRSTIN and NFILE are not consistent. The program will print the 
message RESTART REQUESTED AX STATION XX BUT STORED INFORMATION AT SEQUENCE XX 
IS AT STATION XX and the run will terminate with STOP RESTAR. The user should 
correct the value of IRSTIN or NFILE and rerun the case. 
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Subroutine 

ADI 

ADICUP 

BCONB 

BCONTs’ 

COEFT 

COEFTl 

CON 

COPY 

CROSEC 


DMATl 


DESCRIPTION OF SUBROUTINES 

Name Description 

ADI is a general purpose scalar ADI solver. 

ADICUP is a control routine which specifies the 
calling sequence of routines which set up the 
difference equations for the streamwise momentum, 
thermal energy, turbulence kinetic energy and 
turbulence dissipation equations, specify boundary 
conditions and invert the coefficient matrices. 

BCONB controls and manipulates boundary conditions 
placing them in the desired computational form. 

BCONW performs the actual specification of the 
boundary condition for the streamwise momentum, 
energy and turbulence model equations. 

COEFT sets up the finite difference approximations 
for the streamwise momentum and energy equations 
at interior grid points. 

COEFTl performs a function similar to COEFT but 
for the turbulence model equations. 

CON creates printed contour plots of flow variables, 
if called for. 

COPY reads (writes) the required information for 
a restart from (to) a restart file. 

CROSEC performs numerous functions associated with 
a transverse computational plane. The IROTE option 
over writes the n level solution with the n+1 
level solution prior to computation of the next 
axial step. The lOPTl option is used to compute 
area integrals. The I0PT2 option updates the 
turbulent viscosity. The I0PT3 option sets up 
the secondary flow stream function and vortlcity 
equations. The 10PT4 option writes plot files. 

The 10PT3 option computes the Imposed pressure 
gradients and updates the secondary flow velocities. 
The 10PT6 option sets up the scalar potential 
equation. 

DMATl interfaces coefficient generation subroutines 
with the block matrix inverter, MGAUSS. 
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Subroutine Name 


Description 


DYF 


DYZA 

GAUSS 

GEOTRB 

HINVRS 

IFACE 

INPUTS 

LAPLAC 


DYF computes linearized difference coefficients 
of terms composed of the derivative of the product 
of a known function and the derivative of a flow 
variable. 

DYZA computes linearized difference coefficients of 
a transverse derivative of a flow variable. 

GAUSS is a scalar tridiagonal matrix solver. 

GEOTRB reads the ADD code geometry file and loads 
the FG array with metric Information. 

HINVRS is a matrix inversion routine used by BCONB. 

IFACE Interfaces forward marching subroutines with 
the axlsymmetric potential flow solver. 

INPUTS reads the input for a case and performs 
preliminary set up. 

LAPLAC performs the axlsymmetric potential flow 
solution. 


LENGTH 


MAIN 


MGAUSS 


LENGTH initializes the turbulence quantities 
and updates the turbulent viscosity, as the 
calculation proceeds, given values of k and e. 

MAIN is the main program. It controls the case 
set up, the iteration on the mean viscous pressure 
drop and output. 

MGAUSS is a block tridiagonal matrix solver. 


MINVRS 


HINVRS computes the Inverse of a square matrix. 


MMULT HMULT performs matrix multiplication. 

MSUBT MSUBT performs matrix subtraction. 


MTPROF MTPROF sets up initial profiles from data or using 

automated procedures. 

NTRANA NTRANA and NTRANB transfer dependent variables to 

NTRANB and from mass stroage devices and reorders in->core 

storage of the dependent variables as required by 
the solution procedure. 

OUTPUT OUTPUT controls the printing of flow field variables 

in transverse planes. 

PRINTl FRINTl performs the actual printing of flow field 

variables. 
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Subroutine 

PROFT 

READZ 

RESTAR 

ROBTS 

SECANT 

SECFLO 

SHEAR 

STATUS 

TAU 

TPLOT 

Vise 

WALFUN 


Name Description 

PROFT sets up the boundary layer scaling when using 
the automated starting routines. 

READZ reads dependent variables from mass storage 
devices and loads them into two-dimensional arrays 
for printing. 

RESTAR controls the reading and writing of restart 
files. 

ROBTS sets up the transverse computational grid. 

SECANT is a general secant iteration scheme. 

SECFLO controls the solution of the secondary flow 
equations . 

SHEAR computes the friction velocity compatible 
with the log law velocity profile used in the wall 
function boundary conditions. 

STATUS returns control to the program after an 
NTRAN read or write providing there is not an NTRAN 
error. If an NTRAN error is detected, an error 
termination is performed. 

TAU computes the velocity gradients appearing in 
the turbulence model production terms. 

TPLOT creates printer plots of the axial variation 
of specific variables at a particular transverse 
grid location. 

Vise computes the laminar viscosity. Currently only 
the constant viscosity option is available. 

WALFUN is an auxiliary routine called from BCONW 
when wall function boundary conditions are required. 
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PRIHClPAL FORTRAN SYMBOLS 


NPOINT 

MexlauB number of grid points In y direction. 

NZ 

Maximum number of grid points in z direction, 
must be .LE. NPOINT. 

NX 

Maximum number of streamwise computational steps 
or number of potential lines In geometry 
file. 

MZVAR 

The number of different variables stored In the 
z array. 

NIN 

The maximum number of lines of data for a transverse 
plane vhich may be In core at one time. 

MLEVEL 

The number of levels of the dependent variables 
which are stored in the z array. 

NGEOMV 

The nximber of different geometric variables stored 
In the FG array at each y direction grid point. 

CMUT 

The constant Cp In the' turbulent viscosity 
relationship . 

CCONl 

The constant* C]^ in the turbulence dissipation 
equation. 

CC0N2 

The constant C2 In the turbulence dissipation 
equation. 

SIGK 

The constant Og in the turbulence kinetic energy 
equation. 

SIGE 

The constant o^ in the turbulence dissipation 
equation. 

I2EQMD 

Flag Indicating use of the two^equation turbulence 
model. 

MUT 

z array Index for turbulent viscosity. 

KU 

X array Index for laminar viscosity. 

NTAU 

s array index for the velocity gradient term 


appearing In the turbulence equations production 
ten. 
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NLEN 

NRHO 

NVl 

NV2 

NV3 

NENG 

NTKE 

NDIS 

NVOR 

NSTR 

NUI 

NVI 

NWI 

NCPI 

NPHI 

mi 

NH2 

NH3 

NPI 

NH12 

NH21 

NH31 

NH32 

Y(NPOINT,2) 

X(NX) 


z array index for local length scale, 
z array index for density. 

z array index for streamwise component of velocity, 
z array index for viscous v velocity (=0.). 

z array index for viscous w velocity («0.). 

z array index for total enthalpy difference, 
z array index for turhulnece kinetic energy, 
z array index for turbulence dissipation, 
z array index for streamwise vorticity. 
z array index for secondary flow stream function, 
z array index for inviscid u velocity 
z array index for secondary flow velocity, v^. 

z array index for secondary flow velocity, w^. 

z array index for imposed pressure coefficient, 
z array index for scalar potential. 

FG array index for h-j^. 

FG array index for h 2 * 

FG array index for h 3 - 

FG array index for inviscid pressure coefficient. 

FG array index for Shj/3y. 

FG array index for 3h2/3x. 

FG array index for Sh^/^x. 

FG array index for Sh^/Sy. 

An array containing the distance between transverse' 
grid points in physical space. 

An array containing the computational axial location 
of transverse planes. 
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AP 

Ratio of successive axial spacing for geometric 
progression axial grid. 

YZERO 

Reference length. 

RSEF 

Reference radius, equal to YZERO. 

UE(NX) 

Value of u at NE(l)/2.NE(2)/2 stored for axial plot. 

VE(NX) 

Value of Vg at NE(l)/2,NE(2)/2 stored for axial plot. 

R£(NX) 

Value of p at NE(l)/2,NE(2)/2 stored for axial plot. 

WE (NX) 

Value of Wg at N£(l)/2,NE(2)/2 stored for axial plot. 

JX 

Index for axial station. 

NE(2) 

An array containing the number of grid points in 
each transverse direction. 

NS 

The total number of streamwlse steps to be marched. 

UZERO 

Reference velocity. 

RZERO 

Reference density. 

TZERO 

Reference static temperature. 

XENTR 

Distance from ADD code initial plane to mixer code 
initial plane. 

ILAM 

Laminar turbulent flow flag. 

ICDIF 

A preset vail function boundary condition 
difference weight flag. 

ILAP 

A potential flow solver flag. 

LX 

The number of axial stations in the potential flow 
solution. 

LY 

The number of transverse grid points in the potential 
flow solution. 

LSEC 

Secondary flow option flag. 

MLOBE 

The number of points input to describe the lobe shape 
and boundary layers. 

YLOBE(I.J) 

An array containing the lobe shape and boundary 
layer data points. 
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BLHUB(I,J) 


IHUB 

ISHR 

ILOBE 

IWAKE 

A(MREW,MCOL ,NPOINT) 

Z (MZVAR ,MLEVEL , NPOINT 

AG(NPOINT,9,2) 

AG 7 
AGS 

YS(2,2) 

B 

AB 

KMl 

IJ 

PRL 

PRT 

ABENG 

EZERO 

GAMCON 

CMACH 


An array containing the hub and shroud boundary 
layer data points. 

The number of hub boundary layer data points. 

The number of shroud boundary layer data points. 

A flag used to indicate whether lobe boundary layers 
are to be considered. 

A flag controlling Iteration for a wake parameter in 
the initial boundary layer profiles. 

An array containing the block matrix coefficients 
of the difference equations. 

NIN) An array containing the in-core values of 
dependent variables. 

An array containing the transverse direction 
difference weights. 

l./Ax. 


-l./Ax. 


An array containing the transverse grid limits. 

A coefficient in the dimensionless form of the 
gas law. 

A coefficient in the dimensionless form of the 
gas law. 

The inverse of the Reynolds number. 

A transverse grid index. 

The laminar Prandtl number. 

The turbulent Prandtl number. 

A coefficient in the energy-temperature relationship. 
Reference total energy. 

The quantity y/y-l, where y = the ratio of specific 
heats. 

A reference Mach number. 
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ZE(KR0V,MR0y,2) ' 

ZF(MR0W,MR0W,2) 

ZG(HR0W,MROW,2) 

ZH(HR0W,MR0Wp2) j 

MHIGH 

EPS (2) 

IPLOT 

RG 

CP 

VISCOS 

TSTAG 

GAMMA 

DWT(MROW,4) 

DW 

NCS 

NBCON (MROWp 4 pNSETS) 
NBCONC (MROWp 4 , NSETS ) 
IWALLF 
XCBC 

ITCNT 

ITMAX 

RUFLUX 

FG (NGEOMV , 2 , NPOINT) 
YSAVE(NPOINT,2) 

PEXIT 


Dunmy arrays used in setting up boundary conditions 
for coupled systems of equations. 

Maximum of ME(1) or NE(2) depending on ADI sweep. 
Grid stretching parameters. 

Plot flag. 

Gas constant. 

Specific heat. 

Reference kinematic viscosity. 

Reference stagnation temperature. 

Ratio of specific heats. 

Boundary condition coefficient array* 

Boundary condition source term. 

Variable indicator in boundary condition routines. 
Boundary condition specification array. 

Boundary condition specification array. 

Wall function flag. 

Axial location where a change in boundary conditions 
occurs. 

Viscous pressure drop Iteration index. 

Maximum number of pressure'^drop iterations allowed. 
Reference mass flux. 

An array containing the needed metric Information 
at the current axial location. 

An array containing the computational locations 
of the transverse grid points. 

Pressure used in the thrust calculation. 
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ZVAR(100,4) 

IRSTIN 

IRSTOT 

JRSTIN 

JRSTOT 

NFILE 

NSAVED 

COOR(NPOINT,4) 

COORA(NPOINT) 

USTAR(NPOINT,4) 

VKC 

PCORR(NPOINT) 

IPCOR 

lAXI 

DELMIX 

YCL(NPOINT) 

USTRMl 

USTR^12 

ESTRMl 

ESTRM2 

TUI 

TU2 

ALEN 


Dummy array used in axial plots. 

Restart flag. 

Restart dump flag. 

Restart file unit number. 

Restart dump file unit number. 

Sequence number of restart data. 

The number of restarts saved on file JRSTOT. 

An array containing the physical locations of the 
grid points. 

An array containing the physical value of the axial 
coordinate at the previous axial step. 

An array containing the dimensionless friction 
velocity. 

von Karman's constant. 

A dummy array. 

Potential flow normalization index. 

Geometry option flag. 

Shear layer thickness. 

Dummy array used in setting up lobe shape. 

Dimensionless velocity of the primary (turbine) 
stream. 

Dimensionless velocity of the secondary (fan) 
stream. 

Total energy of the primary stream. 

Total energy of the secondary stream. 

Turbulence intensity of the primary stream. 
Turbulence Intensity of the secondary stream. 

Free stream length scale. 
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USCALE 


Metric scale factor 


ICOMP 

UCOMP(NPOINT) 

IDATA 

DELHUB 

DELSHR 

ICOAX 

IKECUP 

IDIS 

IPRNC20) 

ICON 

DISTL(2,NPOINT,NPOINT) 
ATS (NATS, NPOINT) 


Run node flag. 

A dunny array. 

Starting option flag. 

Hub shear layer thickness. 

Shroud shear layer thickness. 

Starting option flag. 

Turbulence model coupling option flag. 

Turbulence starting option. 

A 20 bite binary output word. 

Printer contour plot option. 

An array containing the distance to a surface 
and the appropriate boundary layer thickness, 
used in starting. 

A dumny array. 
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CODE ADAPTABILITY 


The current MIXER code was written in a form which readily provides a means 
for adapting the code to changing needs. For example, the interface 
with the ADD code is performed in subroutine IFACE for the axisymmetric potential 
solver and in GEOTRB for the basic MIXER calculation. Thus other sources of 
coordinate generation could be used with the MIXER code through proper modifica- 
tion of these subroutines. Similarly, subroutine MTPROF controls the means of 
setting initial conditions and can be modified to provide additional starting 
options and capabilities. 

With regard to adapting the MIXER code to a machine other than a UNIVAC 1110, 
it is noted that the systems routine NTRAN is used to perform most of the reading 
and writing on mass storage. If the code is to be used on a machine which does 
not have this system routine, either the routine must be replaced or a subroutine 
which mimics NTRAN can be written using available system resources. Additionally, 
PARAMETER and INCLUDE statements are used in many areas of the code to specify 
dimensions and to access commons. If these statements cannot be used, they must 
be changed to suit the contemplated machine. 
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Figure Bl, - Flow chart of the Mixer code calculation procedure. 
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Figure Bl.- Continued 
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Figure B3t~ Definition of input parameters for the ICOAX option 
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Figure B4.- Specification of the lobe shape and hub and casing 
boundary layers. 
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Figure B7. - Concluded 
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Figure B9. - Output table of grid and difference weights. 



r 


• • POl^NTl^L FLi>/» SOUiTjoiJ •• 


• •• 

40 j COf*'7F9iiP'5 IW 

1'1 f fF>*M lU'iS, 

FHS.HA I » . 

1 <111110-02 ,1 

3 46 1 H-i| 1 

. mooo-n? 

, M1604.no 





1 

/ 

1 

h 

S 

t-I 

7 

6 

9 

10 


. onnnn 

. 1 

,21i>/0-ii 1 

, N'-ii 1 

.*6 4*0-0 1 

.5006H-I) 1 

,(j46h5-01 

.Hi 737-ni 

.9.1712-01 

.10564*00 

t » 

. niinno 

. 115P*i' i| 

,2 111 r<-ii 

, 1«»'7M-|ll 

.4^ l*|)-il I 

.5<«066-0 1 

,f.466 3-0l 

,61715-01 

.91710-01 

,10564*00 

16 

.0000(1 

. ll‘i.>6-'U 

.2 3fniO-oI . 

, J**r<Hn-i)l 

. ,*(<.1*0-01 . 

. ,5*1066-11 1.. 

.6si8b9-n 1 

.61724-01 

-.93702-01 

.10563*00 

17 

.ooono 

. ns>7-n) 

.2 

,34^h2-01 

,*6 3* J -0) 

.5-106 1-0] 

,69853-01 

.81719-01 

,93640-01 

.10582*0(1 

16 

.00000 

, 1 s?a-ni 

.21064-01 

.1<,M14-^)| 

,*♦, i*i-o) 

.56060-01 

,64645-01 

,61707-01 

,93675-01 

,10560*00 


.00000 

, 1 1529-01 

,23067-01 

, 346H.7-01 

,*#> j*2-(l ) 

.5-105 7-01 

.64636-01 

,61642-01 

.93656-01 

.10576*00 

14 

.00000 

. 1 1510-ni 

.23(W(1-0J 

. 1*690-0 i 

,46 14J-01 

,54053-01 

.64626-01 

.81676-01 

.9.3614-01 

.10575*00 

1 3 

.000 00 

. 1 15 

.2309 1-01 


, *♦> 1*4-1) 1 

.540*4-0 1 

,64616-01 

.61659-01 

,93610-01 

, 10572*00 

1 ^ 

. 0 0 0 0 II 

. 1 I53*-0| 

.2 1046-01 

, 1.^^97-(l j 

,*(. 14-,-d 1 

,56045-01 

.64605-01 

.61641-01 

.93584-01 

.10568*00 

1 1 

.nnooii 

.1 51h-it| 

.2 n (lO-ll 1 

, 1*701-01 

,*63*7-0 j 

,56042-01 

.64795-01 

,61623-01 

,93556-01 

.10565*00 

10 

.. .00000 

..i.,53d-0L— 

^23iil4.-0i... 

.....34/OStOI 

... .46J49-IU . 

. .S804Ut:.Ul..- 

-. 64785-01. 

._al6U4-01 

._,93b2fl-(U ._ 

..10561*00 

9 

.noono 

.1 540-ni 

.21106-0 1 

,3*710-01 

. *6352-01 

.56038-01 

,69777-01 

.fllS87-ni 

,93500-01 

•10557*00 

6 

.00000 

.1 S42«01 

.21112-01 

.3*715-01 

.46156-01 

,56o:i7-oi 

.69769-01 

.81570-01 

.93472-01 

.10552*00 

7 

.00000 

. 1 545-01 

.23] 16-01 

,34721-0) 

,46360-01 

.56038-01 

,69763-01 

.61554-01 

.93*46-01 

.10548*00 

6 

.0000(1 

,1 547-01 

,23l2l-(U 

, 1472+-01 

,46.165-01 

.56040-01 

.69756-01 

.81541-01 

.91421-01 

*l0545*00 

5 

.00000 

.1 54>*-nl 

.2 112^-01 

.34732-01 

.46 171-01 

.56042-01 

,64755-01 

, 81529-01 

.93399-01 

,10541*00 

4 

.00000 

. 1 551-01 

.21129-01 

,14736-01 

,46177-OJ 

.56046-01 

,64/53-01 

,81520-01 

.93380-01 

,10538*00 

3 

,00000 

.1 ISS2-0I 

,23n?-()| 

.34743-01 

,46 162-01 

.56049-01 

,6475 l-Ol 

.61514-01 

.93366-01 

.10536*00 

? . 

- .0000.0 

._Lib53-0L_ 

.-m2il34sM- 

— L34J.42-;.ui. 

.4b38brUl..., 

... 5 3058^0.1 . » 6a.Ii 3s Oi. 

— .815 1 Os.0 1 ..93353 - 01.. 

■•10534*00 

1 

.noooo 

.11553-01 

.23135-01 

.34746-01 

.46JH6-0I 

.58053-01 

,69753-01 

.81509-01 

.93354-01 

•I0533*0n 


.. UBl t 

.I{bo«J*oo 

.hBos^oa 

n Hon*o'> 
79S*00 
.il?9o*on 

...iitsst.oo 

.n779*00 

.11771*00 

.n7Bfl*00 

.ll761*ni) 

, 1 1 7S«*no 
. 1 1 754*00 
. Il75?*m 
. 117.51 *Q1_ 


. . 12 

• 1 ioo7*no 
. j 1007*00 

• I.10«h*00 
,.l30B4.*aO. 

.nooi*oo 

.(1'>77*00 

.l3073*00 

.nOAO*OQ 

.1.3062*00 

.l3n55*00 

,i304H*00 

.I3o?6*oo 

.130lfl*00 

*13011*00 

.13005*00 

.11000*00 

.12096*00 

Lzauij* Qa. 


^ 13 U . 

« 1439 1*00 .15744*00 


15 

. 17 14fi*nn 


.U39?*00 .15743*00 .17145*00 

.14391*00 .|574?*00 ,17143*00 

__Li.4.iHaiouu .4-512 Q:».dq .ii 139*00 

*143B5*00 ,15734*00 .T7133*00 


*143B5*00 
•J43flU*00 
.14374*00 
fel436n*00 
.14360*00 
.4352*00 
, 4343*00 
.^L. oil* OIL 
* 4323*00 
. 4313*00 
•14303*00 
4 14294*00 
. 4205*00 
.14270*00 
.14273*00 


,15734*00 .17133*00 
.157?«*00 .17 25*00 
.15I?0*00 .17 17*00 


. 15720*00 
.15712*00 
.15702*00 
.15692*00 
,15660*00 


16 IT 

.16603*00 .20116*00 

,1H602*00 ,20116*00 

,16599*00 .20113*00 

.,»LUh94 too .2010 7 VOO. 

,1«587*00 ,20099*00 


. I 7117*00 
.17106*00 
.17094*00 
.17061*00 
.17067*00 


*l56b^a0..„-.Uo52*Q0 
,1565^00 .17036*00 
.15642*00 .17020*00 
.15630*00 .17004*00 
.15617*00 .16969*00 


.16576*00 .20069*00 

.16566*00 .20077*00 

.16555*00 .20063*00 

.16541*00 .20047*00 

.16526*00 ,20029*00 

.1H509*00 ,20009*00 

„. 1 54 9.0 ♦ 00 _m.i99.aH* 0 0 . 
.16471*00 .19965*00 


•15630*00 

.15617*00 

,15606*00 

.15597*00 

.15591*00 


.17004*00 

.16969*00 

.16975*00 

,i5963*no 

.16955*00 


.♦ilU — mlb5H9t.0Q m.I69b2t00 .JLa36.l±ilQ ^19038*^0— 


.16471*00 .19965*00 

.16452*00 .19942*00 

.16412*00 .19919*00 

.16413*00 ,19696*00 

.16396*00 .19675*00 

.171.361*00 .19656*00 

.16171*00 ,14843*00 

.rnJLa36.l±im .19038*80- 


1ft 

,21691 *00 
.21690*00 
.21666*00 
...21680.* 00 
.21671*00 
.21660*00 
.21647*00 
.21631 *00 

.2j611*00 

.21544*00 

.21572*00 

-.81548*00 

.21523*00 

.21496*00 

.21469*00 

•21441*00 

.21415*00 

.21392*00 

.21375*00 

. 813694.00 


19 

.23324*00 
.23322*00 
.23318*00 
— ,23348*00. 
.23303*00 
.23292*00 
.23278*00 
.23262*00 
.23244*00 
. 23223*00 
.23200*00 

.—.231.76*00 . 
. 23148*00 
.23119*00 
.123069*00 
,?305B*00 
.23027*00 
.22996*00 

.22975*00 

_*88987.*.00- 


20 

.25012*00 
.25011*00 
.25007*00 
■-^25001*00 
,24993*00 
.24962*00 
•24964*00 
.24954*00 
. 24936*00 
.24917*00 
.24695*0(1 
...24871*00 
. 24645*00 
.24616*00 
.24766*00 
.24757*00 
.24725*00 
.24694*00 
.24667*00 
.-.24658*00 
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Figure Bll. - Sample printout of metric Information. 
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